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DETERMINATION OF METEOROID ENVIRONMENTS FROM 
PHOTOGRAPHIC METEOR DATA 


SUMMARY 

A mathematical model is used to represent 
Opik's 1958 physical theory of meteors in a form 
convenient for programming the computation of 
meteoroid photometric mass values. Sub-samples 
of 333 photographic meteors from McCrosky and 
Posen's sample are selected with respect to 
magnitude scaled for minimum velocity. A statistical 
comparison between the 1958 Opik result and the 
1933 Opik provisional result, the Harvard-Meteor- 
Project basis for mass values since 1938, is given. 
Meteor and orbit parameter distributions and mean 
cumulative flux in absolute units for mass, momen- 
tum, and energy are given separately for the 
terrestrial influx and for the lunar and interplanetary 
vehicle onfluxes. 


INTRODUCTION 

Scope 

Way, Gove, and van Lieshout [1] gave a systems 
functional diagram and description of what might 
appear to be a succinct and sufficient delineation of 
scope category for scientific information. The 
present work will not seem to fit any of their four 
stages: reports or papers, compilations of data or 
tabulations of functions, review articles, and 
treatises. Some aspects of each of those categories 
will be evident. The work does not aspire to the 
level of a treatise in the sense of Opik's [2] or 
Levin's [3] , since it is not so definitive as either; 
e. g. , the physical theory of meteors is treated 
herein mostly indirectly by statistical plausibility 
arguments. Considerable literature is surveyed, 
but with only limited objectives as they bear on a 
vied point or on some incidental product of the 
analysis. 

Any descriptions, explanations, derivations, 
etc. , which are already available on the subject 
are not intended except insofar as they are 
warranted by the particular issue. Textbooks are 
available which approach the problem from a variety 
of special interests; e. g. , for astronomy by Lovell 
[4] , for engineering by McKinley [5] , and for space 


science by Hawkins [6]. Broad general surveys of 
the subject are available which are relatively less 
involved in tedium than this one; e.g. , those recent 
ones by Cosby and Lyle [7] and by Vedder [8] , by 
Dalton [9] in 1962, and our 1963 review [10] and 
annotated bibliography [11]. The reader wanting a 
relatively short account of meteors, but more than 
just a descriptive treatment, will appreciate 
Whipple and Hawkins' [12] excellent 45-page paper. 

A very good and even shorter introduction to meteors 
is given in the first part of a paper by Opik [ 13] . 
Numerous works are available which relate to the 
hazard of meteoroid impact in space; e. g. , con- 
tracted studies (as by Frost [14]) , papers in 
engineering journals ( as by Christman and 
McMillan [15], and astronautical symposia papers 
(as by Dalton [16]). But, toward hazard modeling, 
this report goes only so far as to show the particle 
flux and distributions of mass, velocity, momentum, 
and kinetic energy. 

It must be understood that no inferences can be 
made from this analysis concerning the possibility 
of a flux of orbiting meteroids in cislunar space, 
the possibility for which has been studied theoreti- 
cally by several authors, e.g. , Katasev and 
Kulikova [17] . Any such particle flux would be in 
addition to that inferred from any meteor statistics, 
but would not be independent of satellite impact and 
puncture experiments. There is not sufficient 
information about this possibility to know confidently 
whether it may or may not exist. 

Emphasis is placed on the more pertinent aspects 
of the study, rather than on all aspects equally. 

Some of the results are of little interest per se, but 
have been included primarily to broaden the per- 
spective basis for deciding the relative plausibility 
of the alternative models for particle mass. Such 
results are passed over lightly. Some of the other 
results which are of interest per se, may not be of 
much direct interest for systems engineering. For 
example, the terrestrial influx, the lunar influx, 
and the onflux of particles impacting onto a vehicle 
in interplanetary space are of direct interest for 
systems engineering. But the concentration of 
particles in interplanetary space, independently of 
the probability of interaction with any space vehicle 
or other system of interest, is also helpful for 
tying the results in with other independent information, 
e. g. , zodiacal light and solar F-corona. 



Method 

Instead of using weighting functions or a strati- 
fied analysis to reduce the bias from physical 
selectivity with x'espect to velocity, only those 
meteors are selected for further statistical analysis 
which should have had a reasonable chance of having 
been detected even if the same particle would have 
only the minimum infall (or escape) velocity. A 
sub-sample of such sporadic meteors is selected 
for eacli of two different models of luminous efficiency. 
Statistical analyses with respect to several param- 
eters are made with means, correlations, scatter- 
diagrams, cumulative distributions, etc. The 
x'esults ax’e compax’ed cx’itically in x'elation to those 
reported for various efforts by others. Weighting 
functions are used in part of the analysis, fox - pur- 
poses other than to reduce the bias from physical 
selectivity with x'espect to velocity. However, 
methods to reduce their x’ole, without loss of 
effectiveness, are used. For example, a slightly 
varying weighting of the arithmetic celestial 
latitude of the meteor radiant is substituted for the 
stronger weighting of algebraic values to reduce the 
bias from physical or geometrical selectivity due 
to the northern location of the data*colleeting locale. 

The extensive use of point-by-point cumulative 
distribution or relative cumulative distribution plots 
instead of the moi'e conventional histograms of 
frequency distribution or probability density function 
for a variable of interest may seem strange, especially 
to those with backgrounds in engineering rather than 
physics, mathematics, etc. The historic practice 
of histograms was necessary to minimize compu- 
tation and hand plotting before the prevalence of 
fast automatic digital computers and associated 
machine plotters. Most mathematicians and some 
physicists will prefer the cumulative distribution 
plots as used herein because they contain nxox'e than 
just the information necessax’y to construct the 
histogx'ams. The most common histogram, with 
bars of equal width, can be constructed from a 
cumulative distribution plot by erecting ox'dinates 
at the boundaries of equally spaced intervals of the 
abscissa. The slope of the line segment or chord 
joining the points of intex'section of the cumulative 
distribution plot and the adjacent ox'dinates is the 
height of the bar in the histogram. When the 
relative cumulative distribution plot of a continuously 
varying random variable is fitted with a smooth 
curve (usually not by algebx'aic computation but by 
judging the best fit of an adjustable curve marker) , 
the slope of the curve at any abscissa value is the 


corresponding experimental value of the probability 
density function of the vai'iable. 

Mathematical Considerations 


In addition to scatter diagrams, the programmed 
computations for this analysis ax’e mostly for the 
means, standard deviations, and cumulative 
distributions of, and correlations between, several 
meteor parameters of intex'est. Comparable values 
are computed for each of the two partially over- 
lapping selected sub-samples. These results are 
determined not only with uniform weighting, but 
also with other statistical weighting; e. g. , the mean 
velocity for the lunar influx is computed from the 
same data as for the atmospheric influx, but with 
different statistical weighting. 

The data used in these computations are from 
those of McCrosky and Posen f 1 8 J (1961): 2529 
meteor's, including 463 which were designated as 
belonging to showers, with mass and magnitude 
data missing from 26 of the remaining spox'adics. 

This left 2040 sporadic meteors with values for 
absolute photographic magnitude. The latter wex-e 
extrapolated to correspond with the minimum infall 
velocity before selecting the 333 brightest meteors 
as a sub-sample for statistical analysis. The 
extrapolation depends on the formulation for the 
meteor luminous efficiency. But it was intended 
that two alternative models for luminous efficiency 
should be compared. This selection gave two 
partially overlapping sub-samples, each of 333 
meteors, one sub-sample for each of the two models 
of luminous efficiency. 

Some of the meteors, which McCrosky and 
Posen ( IS] said wex'e due to particles in heliocentx'ic 
ox'bits of lax'ge aphelia, and which they called 
"hypex'bolics, " had no values tabulated for two of 
the twelve parameters considered in this analysis 
(aphelion q 1 and comet-astex’oid critex'ion K) . 
Sub-sample "A" (in consideration of Opik's [19] 

1933 meteor physical theox-y ) includes two 
"hypex’bolics." Sub-sample "B" (for Opik's [2] 

1958 x-esults) includes 11 "hypex'bolics. " There 
ax-e 274 members common to both "A" and "B, " 
including the two "hyperbolics" of "A." 

The error's which McCrosky and Posen [18] 
x'eported for the elements of the meteors by the 
graphical reduction (e.g. , three percent rms 
deviation of velocity) , and the errors in the reported 



values of absolute photographic magnitude M 

Po 

according to a subsequent study by Kresak [20J 
(0. 25 probable error) , are small enough to be 
ignored in this study. The uncertaintaies which 
play the main role in the propagation of confidence 
in the indicated calculations are those due to the 
finite size of the sample. The size of the sub-samples 
(333) is a compromise intended to foreshorten any 
physical selection effects and to avoid the compli- 
cations of a stratified analysis. 


For an analysis of two parameters x and y in 
a sample of N members, with parameter and weight- 
ing values x., y., and w., with w. adjusted so that 


N 



i=l 


( 1 ) 


the mean x and standard deviation <j^, for the 

parameter x, and the correlation r between x 
, xy 

and y are 


N 


x = (1/N) £ 
i=i 


w. x. 
i l 


a 

x 


N 

(i/N) y w.(x. 

Si 1 1 



(2) 


(3) 


r 

xy 


N 

(-' cr a ) ~ Yj w (x - x) (y - y) , 


(4) 


where, by equation (1), the weighting values w 

have been normalized so that N is both the number 
of elements and the total statistical weight. Here 
the number of elements N is 333, except in calcula- 
tions involving either q' or K, where it is 333-2 and 
333-11 in sub-samples "A" and "3, " respectively. 
Equations (1) through (4) differ from textbook 
formulations (e.g. , Reference 21) only for w. £ 1. 

The confidence intervals for the means and corre- 
lations are calculated in the usual manner, i. e. , 
x is approximately normally distributed with 
standard deviation u^/n/N, and the parameter 

2 log |(1 + r X y)/ (1 - i’ ) ] Is approximately 

normally distributed with standard deviation 

_ i 

(N - 3) For example, when N is 333 and the 


computed correlation is zero, the probable error of 
the correlation is 0. 037. However, the textbook 
formulation for the confidence intervals in r , 

xy 


although unavoidable if it were necessary to 

establish percentiles of very high order, can be 

further approximated for the present purposes. 

These results were checked for a matrix of 

representative values of N and r over the 

xy 

intervals 59 s N < 333 and -1 s r £ 1, and they 

xy 


have shown that r is approximately normally 

distributed throughout the interquartile range, with 
standard deviation 



2 (1 - | r I ) (N - 3)"* 


for |r | S: 0. 9 
xy 

( nv\ °- 87 
l cos — J (N - 3) ' 

for |r 1 < 0. 9 
xy 


(5) 


For example, the semi-interquartile range (or 
probable error) of r is found, by multiplying 

equation (5) by 0. 6745, to be accurate to within 
0.001 for |r 1 £ 0. 3 and to within 0.002 for 


xy 

|r | > 0.3, in comparison with the more tedious 
xy 

transformation. 


Confidence intervals for the computed values 
of the cumulative distributions of the parameters of 
interest can be approximated by a formulation 
previously given by Dalton [22] for uniformly 
weighted data. The normal approximation, though 
expedient, is grossly inaccurate. The binomial 
formulation for the percentile p^ of order C 

(0^C<1) for the relative cumulative distribution 
(e. g. , Reference 23) is implicit, rather than 
explicit, with the result that users' results are 
usually read from tabulations (e.g. , reference 24) , 
subject to the errors of interpolation. When the 
sub-sample "A" or "B" has been ranked with respect 
to a parameter of interest, and a value of that 
parameter has been chosen to divide the sample 
so that nj. s N/2 of the members are on a side of 

it, then the confidence C (1/2 :£C < 11 that the 

fraction of the population on that side does not 

exceed p is related to o by 
s * s 


3 



P s = 1 

whe re 


(1 - C)/e 


n f 

S w > 

1=0 


-1 [(3+ 20) M] 


' " f V 

N - ) w /2 

^ ' J 


( 6 ) 


f 


6 = 0 when V = 0 
i=0 


f 

1 when Yj # 0 . 
i=0 


For example, with N = 333 and n = 0, the values 

1 1 3 

of p corresponding to C = — and C =— are 0.0021 
S 2 4 

"f 

and 0.0042, respectively; or if ^ w. were 33.3, 

i=l 1 

the median and third quartiie of p^ would be 0. 10? 

-. 0. 122, respectively. 

When calculated values for sub-sample "A" are 

being compared with the corresponding results for 

"B, " it is necessary to account for the 274 members 

common to both ’’A" and "B, " and for their role in 

equations (1) through (6), before estimating the 

statistical significance of the differences between 

results. The variances of the parameters were not 

computed separately for the N 0 members which are 

common to both sub-samples "A" and "B. " But in 

computing the confidence for the difference between 

corresponding values from the two sub-samples, the 

mean variance for the N. members of "A" not in 
A 

"B, " and lor the N members of "B" not in "A, " 

B 

can be approximated by the mean variance for the 
corresponding total sub-sample. Then the 
difference between the mean x ^ of a parameter 

x in sub-sample "B" and the corresponding mean 

x A front sub-sample "A" is approximately normally 

distributed, with standard deviation 



•v 


= (1/Nl 




(71 


For example, with uniform weighting, the values 
for the mean air-entry velocity for "A" and "B" are 
16. 62 and 20.23 knt/sec, respectively; and the 


corresponding values for the standard deviation of 

the air-entry velocity are 5. 01 and 10.32 km/sec, 

respectively. In this case, since the data for the 

"hyperbolics" also enter the computation, N is 

333 and both N. and il are 59. Then the calcu- 
A B 

lated value and standard deviation of the difference 
in uniformly weighted mean velocity are 3. 61 and 
0.27 km/sec. Similarly, in terms of the terrestrial 
statistical weighting and probable error, the mean 
air-entry velocity for "B" is 18. 83 t 0. 31 km/sec, 
which is higher than that for "A” by 2. 66 i 0. 15 
km/sec. 

With the same approximations which were used 
in deriving equation (7) , the difference between the 
"A"- and "B''-calculations for a particular corre- 
lation can be expressed by 


r - r « ( 1/ Nj 
Xy B xy A 


N fr ) 

B xy N 


- N ( r ) . . 
A 1 xy' N 


( 8 ) 


where (r ) is distinguished from r 
the latter is the value of r in "B" with all 


in that 
members 


xy 

N, whereas the former is that in "B" with only 
members N and similarly for the subscript "A. " 


Since, by equation (5) , equation (8) is a linear 
function of approximately normally distributed 
variables, it is also approximately normally distrib- 
uted. with standard deviation 


,[-»(' '*(' NJ)‘‘ 

* 

] 

->J 4 J 


r r,J / 



^i(i >') 

1 

4 «A-» 

. 

° A 7 * 1 
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For example, with terrestrial statistical weighting, 
the correlation between air-entry velocity and 
"Model-A" mass m^ in sub-sample "A” is 0. 020 

± 0. 037 while the correlation between air-entry 

velocity and ”Model-B" mass m in sub-sample 

B 

"B" is 0. 191 ± 0. 035. In consideration of the over- 
lap between "A" and "B, " the difference between 


4 



these correlations is 0. 171, and to the extent that 
one might neglect the fact that "Model-A" mass and 
"Model- B" mass are not really the same parameter 
and may play a somewhat different role in the common 
members N 0 , the probable error of the difference 
between these correlations is 0. 009. In other words, 
the "A" and "B" values of the correlation between 
mass and velocity differ quite significantly from a 
statistical point of view. 

Not only is the same analysis for sub-sample 

"A" repeated for "B, " but also, in addition to 

uniform weighting, the analysis is repeated three 

more times, each time with a different function for 

the statistical weighting w: terrestrial weighting 

w , lunar weighting w and spatial weighting w . 
t 1 s 

Just as the two sub-samples, "A" and "B, " are 
not statistically independent, but, having some 
common members, must be interpreted with 
equations (7) through (9) , so also do any two of 
the analyses with different functions for the statisti- 
cal weighting w. The following summations play 
essentially the same role as N^ and 


993 

) |w - 11 = CO. 5 for sub-sample "A" 
i=i l 

= 09. 3 for sub-sample "B" 


333 

Y, lw - w I = 120.3 
i=l i i 

= 117. 9 
333 

V )w - w, I = 193. 5 
h 1 . s. 1. 

i=l l l 

= 206. 9 


for sub-sample "A" 
for sub-sample "B" 

for sub-sample "A" 
for sub-sample "B" 


For example, with terrestrial weighting w , the 

mean and standard deviation of the geocentric 
velocities of the particles are 14.47 and 9.88 km/sec, 
respectively, for sub-sample "B, " while with lunar 
weighting, w^, they are 18.04 and 11.09 km/ sec, 

respectively. These means could be reported as 

14. 47 ± 0. 37 and 18. 04 ± 0. 41 km/ sec; and, by 

substituting 117. 9 for N and N in equation (7), 

A B 

the difference between these statistically dependent 
mean values would be 3. 57 ± 0. 33. The fact that 
these means for the geocentric velocities in the 
terrestrial and lunar influxes show any statistically 
significant differences is due to the velocity 


dependence of the gravitationally modified particle 
influx cross section of the earth. 


The particle influx cross section which the 
earth would have without a gravitational field, 
relative to that which it actually has, is proportional 
to the square of the ratio of the geocentric and air- 
entry velocity, with a different value for each parti- 
cle. But when this result is multiplied by the 
terrestrial weighting w^ and summed over all 333 
particles, one gets the number of particles, within 
the same mass interval, which should have 
encountered the earth without gravitational focusing. 
When this result is divided by 333, the size of "A" 
or "B, " one gets a good approximation to the mass- 
cumulative influx of the moon k^ (per unit area) 

relative to that of the earth (within the mass interval 
of the data) ; i. e. , 



= 0.470 for sub-sample "A" 


( 10 ) 


= 0. 540 for sub-sample "B" 
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DISCUSSION 

Meteor Luminous Efficiency 

A meteor particle (or meteoroid) enters the 
atmosphere with some initial (or air-entry) velocity, 
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V , essentially corresponding with the meteor 

beginning height H . Some fraction of the air-entry 

kinetic energy is converted to meteor radiative energy. 
Following Opik [19]. who calculated this entity in 
1933, the integrated meteor radiative energy in the 
wavelength interval 4500 to 5700 A, relative to the 
air-entry kinetic energy, is called the meteor 
luminous efficiency p. t)pik [19] gave values for 
two components of p, a component j3 0 for impact 
radiation, and a component p^ for temperature 

radiation. He concluded, "In most of the meteors 
visible to the naked eye, temperature radiation is 
insignificant compared with impact radiation. " A 
plot of his tables of values for jin versus V is 

oC 

shown in Figure 1, separately for visual and 
telescopic (i.e. , fainter) meteors. He suggested 
that the transition between these (quite different 
visual versus telescopic trends) might begin at the 
fainter visual meteors. In 193S, Whipple [25] 
used those results as the basis of his model for 
visual meteors, "Model-A" in this report: 



(H) 

t 0v = S. 5 x 10 -a sec/km , ^ 

where the constant t<j is usually called the 

luminosity coefficient (e.g. , by McCrosky and 
Posen 1 18]), visual, in contradistinction to to , 

photographic (e.g., by Verniani [26]). This 
"Model-A" was also used by Whipple [27] in 1943 
for his meteor studies relating to the properties of 
the upper atmosphere. Indeed, as Whipple [27] 
stated, the results reported in 1933 by Opik [19] 
are well represented by "Model-A. " But it is at 
considerable disparity to the more studied analysis 
published by Opik ]2] in 1958 and illustrated by 
tables of calculated values also by Opik [2S] in 1963. 

When the product of the meteor luminous 
efficiency (visual or photographic) and air-entry 
kinetic energy' is equated to the integrated intensity 
(visual or photographic) , a mass value can be 
calculated which is called photometric mass. Such 
values were tabulated for the Super-Schmidt photo- 
graphic meteors by McCrosky and Posen [18]. 

They used the "Model-A" luminous efficiency and 
derived the integrated intensities by Hawkins' [29] 
short-trail method. It must be understood that, in 


order for them to use "Model-A, " they converted 
the photographic intensity to an equivalent visual 
intensity by means of a color index, as discussed by 
Verniani [26] (i.e. , the relation between radiative 
energy and meteor intensity is somewhat different 
for visual intensity and photographic intensity). 
Another value for the mass of a meteoroid, called 
the dynamic mass, can be calculated by considering 
the measured deceleration of the particle during its 
interaction with the atmosphere. But Opik [2] says 
that the photometric mass values are more reliable; 
and Verniani [26] states that all Harvard meteoric 
masses have been computed with the "Model-A" 
luminous efficiency (i.e. , also photometric). 

While claiming no high accuracy for "Model-A, " 
some authors (e.g., McCrosky and Posen (IS]) 
imply that the tabulated values can be corrected 
when the appropriate scaling factor becomes better 
established; i.e. , as illustrated by Clough and 
Lieblein ( 30] , the scaling factor is the ratio of the 
old and new values of the luminosity coefficient, 

7 0 inequation (11). This impression probably is 

encouraged unwittingly by the practice of expressing 
particle mass as twice the product of the integrated 
intensity divided bv ?« V 3 , where r n V' 3 might 

V V SO V V so 

be better written as 8 , V 2 . Then, when an alternative 
A 

"Model-B" gives a luminous efficiency p the 

sealing factor is p //} which one should not 
A B 

expect to be independent of either mass or velocity. 
Anyway, the product of luminous efficiency and 
photometric mass value is invariant between alterna- 
tive mathematical models: i.e. , 



Thus, what is actually known about any meteor is 
neither the luminous efficiency nor the mass, but 
the product of the two. 

Papers by Kallmann [31] and Opik (32] at the 
1954 symposium in Liege suggested an interpretation 
of Opik's [19] results quite different from Whipple's 
[25, 27] "Model-A." Noting the latter model, 
Kallmann [31] said that it must be assumed that the 
luminous efficiency is not only a function of the 
velocity but is also a function of the size or magni- 
tude of the particle, and that the same must be true 
for the visible path length of the meteoroid. Opik 
[32] noted some data for shower meteors which 
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seem to indicate that the luminous efficiency would 
vary inversely with tire 0.87-power of the velocity. 

Other papers by Kallmann [33] and Opik [34] 
in the 1954 symposium at Jodrell Bank carried the 
disagreement with "Model-A" somewhat further. 
Kallmann [33] used Opik 1 s [19] relation between 
luminous intensity and magnitude to estimate the 
masses of meteors of different visual magnitudes 
greater than zero and of different velocities. He 
showed that the luminous efficiency is a function 
not only of the velocity of the meteor for all magni- 
tudes, but also of the visual magnitude or physical 
size of the particle itself. This was done by showing 
that the duration of visibility of the particle is one 
function of the luminous efficiency and another 
function of the visual magnitude. Opik [34] said 
that from his 1933 work [19] , which he referred 
to as (lie provisional theory, with certain corrections, 
the most probable expression for the luminous 
efficiency was assumed to be 0. 020/V . This value 

would equal that given by "Model-A" at 24 km/sec. 

By 1957, Weiss [35] was led, by examination of the 
literature, to adopt "Model-A" for meteors with 
negative absolute visual magnitudes; but for positive 
magnitudes he chose to equate luminous efficiency 
to 0. 021/V 0 - 3 . 

OO 

A description and technical review of Opik's 
(2) (1958) more definitive "Physics of Meteor 
Flight in the Atmosphere" were given by Eshleman 
[36] . The results remain rather involved and 
tedious to apply; but Opik [28] (1963) mitigated 
that difficulty with tables of values of the luminous 
efficiency (3 and the components /3 0 for impact 
radiation, ( 3 ^ for temperature radiation, and /3^ 

for black-body radiation — each for an array of 
values for the air-entry velocity V and mass m. 

The values for , the only component considered 
toward the construction of "Model-A, " are somewhat 
different from the values which Opik [19] tabulated 
in 1933, and some of the values of /3 (and even of 

/ 3 ^) are relatively quite appreciable. The tables are 

in three parts, separately for stone dustballs, com- 
pact iron, and compact stone pax'ticles. 

Curves for the luminous efficiencies of compact 
iron and of stone meteoroids, relative to the luminous 
efficiency of a stone dustball of the same mass and 
velocity, are plotted from Opik's [28] values in 
Figures 2 and 3. The curves for 7. 4 km/sec are of 
more interest for re-entry of orbited debris than for 


infalling natural meteors. Otherwise, the effective 
similarity between Figures 2 and 3 is quite 
apparent, indicating that nearly the same mass 
would be inferred whether a particular meteor was 
considered to be compact iron or compact stone. 
Figure 3 shows further that, for a meteoroid of any 
mass and any velocity, a higher luminous efficiency 
and consequently lower mass is inferred by the 
assumption of a compact, rather than dustball, 
particle. Therefore, not only the mass-cumulative 
influx, but also, by the liypervelocity impact and 
puncture model which Dalton [37] gave in 1966 and 
presented again after some further refinement in 
1967 [16] , the puncture hazard which one infers 
from meteor data is higher when the particles are 
presumed to be dustballs. 

Results reported by different authors back 
through 1954 are at variance concerning the struc- 
tural integrity and density of the meteoroids 
responsible for the meteor phenomena. Opik [13] 
presents a case for the rarity of compact iron 
versus compact stone meteoroids. Also Opik [34, 

32, 2] reports that the meteors of the visible range, 
as a rule, ax - e not stony compact bodies, but some 
kind of loose aggi’egates of dust particles, "stone- 
flakes" or "dustballs. " Smith [38] finds, by analysis 
of meteor flares, that such dustballs consist of 
thousands of fragments, each about 10 -6 grams. 
Fragmentation, another observational evidence 
supporting the dustball concept, has been reported 
by Jacchia [39] and by McCrosky [40], Meteor 
particle fragmentation is supported on theoretical 
grounds in consort by Whipple's [41] icy comet 
model and by analysis of meteoroid orbits with 
respect to Whipple's [42] comet-asteroid criterion. 
But some of the recent l-esults are contrariwise. 
Ceplecha's [43, 44, 45] analysis of meteor height 
and orbit data indicates three distinct groups — 

A , B and C — with material specific gravities 

v v U 

4.0, 2.2, and 1.4, respectively, with the A group 

and B^ group being most abundant and least 

abundant, l'espectively, in the mass interval from 
about 1 to 10~ 4 gx'ams. Laboi-atory observations of 
the effects of high enthalpy electric arc jet winds 
over exposed rock samples at low ambient pressure 
by Allen and Baldwin [46, 47] suggest that frothing 
and sluffing of compact stone meteoroids may be a 
more general explanation for the anomalously low 
material density. Also, Lebedinets and Portnyagin 
[48] reported a theoretical analysis suggesting that 
the observed pi'operties of faint photographic meteors 
can be interpreted on the basis of a dense- meteoroid 



model, without invoking the hypothesis of an excep- 
tionally friable structure. Therefore, to approxi- 
mate Opik's [2, 28] results, for programming pur- 
poses, a "Model-C" was constructed for compact 
stone meteoroids. This is in addition to the "Model- 
B" for Opik's [2, 28] stone dustballs, a detailed 
derivation for which was published elsewhere by 
Dalton [49] . 

"Model-B" for stone dustballs, which is illus- 
trated in Figure 4 with Opik's [28] values super- 
imposed, relates luminous efficiency p , air-entry 

i3 

mass m_ and velocity V as follows: /3 and ni 
Jd 00 J3 i3 

are functionally independent when either 
V < 14. 8 


and 

log log (mg/3-g) -1 a 0. 0125 + 0. 665 

or 


V a 14. 8 


and 

log log (m /3 ) -1 a - 0. 00304 V + 0. 895; 

13 13 00 

otherwise, 

log log (j 8 b ) _1 = a 0 + a x log log (m B 0 B ) _I , (13) 
where 


a*, = 0. 4042 + 0. 193 (log log - 0. 0682) 
for V a 14. 8 

CO 

= 0. 4042 + 0. 33 (0. 0682 - log log V ) 

CO 

for V < 14. 8 


and 


a, = 0. 0840 + 0. 00113 |V -14.81. 

Correspondingly, "Model-C" for compact stone 
meteoroids is illustrated in Figure 5 with Opik's 
[28] compact stone values superimposed. All of the 
members of sub-sample "B" map into the region of 
functional interdependence of )3 ^ and m^,: 

log log (0 C ) _1 = a 2 + a 3 [log log (m c P c ) -1 

- 0.364] , (14) 

where 


a, = 0. 313 + 0. 079 log V when V > 14. 8 

= 0. 695 - 0. 248 log V when V £ 14. 8 

CO CO 

a 3 = 0.1527 + 0.001315 V . 

The computation of m_ and m_ by equations 

(12) and (13) is facilitated by the fact, as already 
noted, that the product of luminous efficiency p anc 
mass m is invariant from one photometric mass 
model to another, and the values for "Model-A" 
have been tabulated by McCrosky and Posen [18] . 
After is found, it is divided into m^/lg to find 

m , and similarly for p and m . 

Compilation of Representative Statistics 

As recently as 1967, Cook [50] concluded that 
the best research data concerning meteor physical 
theories comes from photographic meteors. The 
photographic Super-Schmidt meteor data of McCrosky 
and Posen [18] give the most extensive photographic 
observational material for 2040 sporadic meteors 
with tabulated values for "Model-A" mass and various 
other data. Their data have been used as the main 
observational material for statistical studies by 
Ceplecha [44], Dohnanyi [51], and Miller [52], 
Ceplecha [44] noted that the limiting and maximum 
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magnitudes of the meteors have a statistical distri- 
bution that indicates a very good homogeneity of the 
McCrosky and Posen [18] material. 

An analysis of the velocity dependence of 
luminous efficiency with Hawkins and Southworth's 
[53] reputedly more accurate data for 285 sporadic 
photographic meteors was reported in 1965-66 by 
Dalton [54, 55]. The results were given with 
different combinations of statistical weighting and 
luminous efficiency functions. The results now 
seem invalid. It is evident now that the sample size 
was too small. Ordinarily one should consider that 
a sample with 285 members would constitute fairly 
reliable statistics, but it is necessary to consider the 
role of the statistical weighting in reducing the 
effective size of the sample. Another uncertainty 
with Dalton's [54, 55] results is due to the role of 
velocity in the weighting function. The velocity 
dependence of the weighting function was adapted 
from Whipple's [42] cosmic weight. Both of these 
disadvantages are avoided in the present analysis 
by choosing sub-samples which do not necessitate 
a statistical weighting as a function of velocity. 


If one is to avoid velocity-dependent statistical 
weighting, then legitimate statistics of meteors can 
be made only when one eliminates the data for those 
meteors which would have been too faint for 
reasonable detectability if they would have had only 
the minimum velocity of interest (e. g. , infall or 
escape velocity, 11 km/sec) . Such a procedure was 
outlined in 1967 by Dalton [49] . It is based on maxi- 
mum effectively detectable absolute photographic 
magnitude M ^ as a criterion, a velocity extrapola- 


tion of the actual absolute photographic magnitude 
M . Jacchia, Verniani, and Briggs [56] found, 


by an analysis of data for 413 precisely reduced 
photographic meteors, that the integrated photo- 
graphic intensity, and therefore also the product of 
luminous efficiency and kinetic energy, is propor- 
tional to ii-il/v 0 - 89 cos 2 / 3 Z The angle Z is 
pm » R R 

the zenith-to-radiant deviation and I is the 

pm 

maximum photographic luminous intensity, related 

to M by 
pm 


M = - (~) log I 
pm V 2 / pm 


Therefore, the criterion for selecting sub-samples 
relatively free of bias due to physical selectivity 


with respect to velocity is 


M = M + 2. 25 log 
Pgo pg 



+ 6. 50 log 



( 15 ) 


where f3 and /Jjj are the luminous efficiency at 
velocity and the luminous efficiency at 11 km/sec, 

respectively. 

The order or rank with respect to M for 

Pgo 

McCrosky and Posen's [18] 2040 sporadic meteors 
depends on whether P/Pu in equation (15) is evalu- 
ated with "Model-A, " equation (11) , or with "Model- 
B, " equation (13). Sub-sample "A" contains the 
333 meteors for which M < 2. 3 by "Model-A," 

pgo 

and sub-sample "B" contains the 333 meteors for 

which M < 2. 0 by "Model-B. " 

Pgo 

The most obvious comparisons between the 
two selected sub-samples are illustrated by the 
scatter-diagrams. Figure 6 shows "Model-A" 
versus "Model-B" mass values, and Figure 7 shows 
mass versus velocity. For example, considerably 
more of the high velocity meteors are represented 
in sub-sample "B, " 


Statistical Weighting 


According to formulas for the transformation 

of coordinates, given by Lovell [4] and Ericke [57], 

a meteor radiant with right ascension a and 

declination 6 has celestial latitude B , where 

e 


sin B =0. 91741 sin o - 0. 39795 cos 6 sin a. 
e 

(16) 

In Figure 8, where the values of (3 for both sub- 
samples are plotted, the values of p are pre- 
dominately positive. 

Kells, Kern, and Bland [58] state that the 
zenith right ascension a ^ of a place is the sidereal 

time at that place; and Russell, Dugan, and Stewart 
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[59] state that sidereal time is usually correct 

within four or five minutes when approximated from 

civil time by assuming that on September 22 the two 

times agree and that sidereal time gains two hours 

each month, and proportionately for each day and 

hour. For the Harvard stations in New Mexico, 

the locale of the meteor data, Hawkins and 

Southworth [53] state that, to obtain the time of 

appearance in local mean tame at the Dona Ana 

station, 0. 29667 must be subtracted from the date 

given in universal time. Also Hawkins [29] gave 

for that station the zenith declination 6 as 

z 

6 z = sin -1 0. 8433343 = cos" 1 0. 5373893. (17) 


With this information, together with equation (16) , 

the zenith celestial latitude B can be formulated. 

z 

Local (Dona Ana station) sidereal day is 


t . — t . + 
ds du 


(t 




365.26 


- 0.29667 


/ 366. 26 \ 
\3G5.26 ) 


(18) 


where t and t, are the numbers of the month 
mu du 

and day universal time, respectively. The zenith 
right ascension cv is given in radians by 


a = 2 7T t . 
z s ' 


(19) 


where t is the local sidereal time expressed as a 
s 

decimal fraction of a day; i. e. , t is the decimal 


residue which is left after subtracting the next lower 
integral value (postive, nil, or negative) from the 
local sidereal day t j , equation (18). By equations 

(16), (17), and (19), the zenith celestial latitude 

B is related to the local sidereal time t by 
p z s 

sin B = 0. 77368 - 0. 213S5 sin 2 

z 

34. 0° K j 3 4 80. 9° 



The values of p versus p for the two sub-samples 
z e 

are plotted as scatter diagrams in Figure 8. 

One should expect the reason the distribution of 
is not symmetric with respect to zero (the 

ecliptic; is related to the fact that the distribution of 
p also is not symmetric with respect to zero. These 

data give no basis for supposing that any more than 
random deviation from a distribution of p symmetric 

about zero should be found from data for which the 
time and locale would have been chosen for the 
ecliptic to be through the zenith. Contrariwise, 
radio meteor data give some evidence for the essen- 
tially symmetric distribution of radiants with respect 
to celestial latitude. Correspondingly weighted 
radio meteor results, from the northern hemisphere 
by Davies and Gill [60] and from the southern 
hemisphere by Keay [61], showed that correction 
for observational selection reduced the earth-way 
apex mode of radiant density well below two accentu- 
ated broad modes at about 60 degrees on either side 
of the apex to the earth's way. Davies and Gill [60] 
found that correcting their radio meteor data for 
observational selection essentially produced symmetry 
in the distribution of radiants in ecliptic latitude, 
and they considered it as evidence for the correctness 
of the statistical weighting. 


10 



I 


The radiant density per square degree of the 
celestial hemisphere is proportional to the probabil- 
ity density function for sin p . 1 The cumulative 

distribution of the positive values of jS , with 

weighting inversely proportional to the square of the 
beginning height H and relative to their total 

.D 

statistical weight, is plotted in Figure 9. The 
weighting is intended to correct for differences in 
the effective area of surveillance of the cameras. 
Also in Figure 9, with the same ordinate scale, is 
the corresponding distribution of the negative values. 
The fact that the number of positive values of p^ 

exceeds the number of negative values by a factor 
of nearly 3 for both sub-samples "A" and "B" was 
already apparent in Figure 8. But Figure 9 suggests 
more clearly how a further weighting function might 
be used to minimize bias with respect to p^. 

The data show that most of the meteors have 
radiants above the horizon. Then, as a tentative 
step, it seems appropriate that meteors which have 
radiants on circles of celestial latitude with not 
more than 90 percent occlusion (by the horizon) 
should be weighted by a factor w in near inverse 

proportion to the apparent (above the horizon) 
fraction. The weighting factor w is related to 

the difference 2X^ between the celestial longitudes 

of the horizon points on the circle of celestial 

latitude B through the meteor radiant. 
e 

The pole of the ecliptic, the zenith, and either 
one of the horizon points, on the circle of celestial 


latitude p^, form a spherical triangle such that the 

angle at the pole of the ecliptic is 1S0° -,\ with 

opposite side 90° , and with adjacent sides 90° -P 

and 90° -13 , where (3 is the celestial latitude of 
z z 

the zenith. The solution to the spherical triangle 
gives 

\ = 0 when fi £ 90' - fi 

e 'e 'z 

cos A. = tan (3 tan 13 when B < 90° -p 
e e z e z 

and tan (3 tan (3 > - 1 
e z 

= - 1 when tan fi tan B < - 1 

e z 

At the Dona Ana station, tan p is always 

positive, 0 < £ 90° when p > 0, 90° < s 180° 

when B <0, and the circles of celestial latitude 
e 

with p > 90° -j3 are entirely above the horizon. 
Then 6 

w = 180°/ (180° -X ) when p < 90° -p 
c e e z 

and 180°/(180°-\ e ) < 10 

= 10 when p < 90° -p 

e z 

and 180°/ ( 180° -A^) > 10 . 

(21) 


1. Let y be some function of p with probability density function f(y) proportional to the radiant density 
per square degree of the celestial sphere, i.e. , 


f(y) 


dN 

dA 


f(y) dy 

cos p dp 
e e 


where dA is the measure of square degrees in the differential element of width dp ^ along the 2« cos /? length 

of the circle of celestial latitude through the radiant, and dh! is the differential number of radiants in dA and 
equals f(y) dy. Therefore, 


dy ~ cos p^ d p^ 

y ~ sin p ■ 

J e 
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The effect of weighting ~ w H 2 instead of 

C B 

~ H“ 2 , is to transform the results of Figure 9 into 

those of Figure 10. The results show that the extra 
weighting from equation (21) is quite effective in 

reducing the disparity between the positive and 
negative branches of the parameter of interest. 
Actually, there is a slight over-correction, which 
could be avoided by choosing a lower limiting value 
for w in equation (21) , say S instead of 10. 

From Figure 10 one gets the impression that 
weighting ~ 1I^ 2 w should be quite appropriate for 

each, member of the sub-samples. But when the 
positive and negative branches of 0 are combined 

with that weighting in Figure 11, the undesirability 
of the result is more obvious. Too much of the 
statistical weight was shifted to too few members, 
leaving a very jagged distribution. But, by looking 
back at Figure 10, one sees that the positive branch 
remains quite uniform even with the weighting 

~ H~ 2 w So, one would accept the positive branch 

of Figure 10 as being sufficiently free of bias with 
respect to But instead of using the weighting 

w , a milder weighting w^. is used to correct the 

bias in 1/3^1 rather than m /i^. By dropping the 

weighting w from Figure 11, and using only ~ H~ 2 , 

one gets the results shown in Figure 12. Then, 
because only l/3 e l and not J is considered in the 

rest of the analysis, the milder weighting w^. which 

can replace w is the ratio of die slopes of the 
c 

corresponding curves in Figures 10 and 12, plotted 
in Figure 13. The values of the ratios for w^. plotted 

in Figure 13 indicate an inappreciable difference 
between sub-samples "A" and "B"; therefore, the 
following mean result can be used for both: 


for l/igl, implicitly involves no weighting with 

respect to the radiant zenith angle Z^. However, 

the weighting ~H“ 2 involves Z because the 
B K 

beginning height H equals the product of the range 

and cos Z^. But McKinley [62] stated in 1951, 

and Opik [32, 63] found also in 1955 by the data 
from the Arizona Expedition, that the meteor 
detection probability function for a visual observer 
is cos Z^, not cos 2 Z^. This means that if the 

observer sees n meteors per hour from a radiant 
overhead, he would see n cos Z^ meteors per hour 

from the same radiant if it were at a zenith angle 
Z R . This suggests that weighting ~ H“ 2 may over- 
correct for Z^, and that such an over-correction 

would have to be adjusted by letting the weighting 

w depend on Z . For example, one would not 
c K 

have been surprised if, instead of calculating the 
length of the apparent segment of the circle of 
celestial latitude, one might have needed to integrate 
some function of along that segment. But the 

weighting in Figure 10 seems to have been sufficiently 
effective in restoring the essential symmetry of the 
radiant distribution with respect to the ecliptic. 

The statistical analysis in the next part of this 
section was done first with uniform weighting. Then 
it was repeated, once with w for the (world-wide) 

terrestrial influx, again with w^ for the lunar influx, 
and with w for the population of particles in inter- 
planetary space regardless of the probability of 
collision with the earth, the moon or an interplanetary 
vehicle. The weighting functions, with coefficients 
adjusted so that the sum of the weights equals the 
sub-sample size 333 in accordance with equation (1), 
are 


1. 10 

for 0 := sin 

'*e 

i < 

0. 32 

0. 61 

for 0. 32 < 

sin 

,P e 

| <0.45 

0. 94 

for 0.45 == 

sin 

‘^e 

1 < 0. SO 

1.22 

for 0. 80 £ 

sin 

W e 

I 1 


( 22 ) 


The equation (21) weighting w for fS the 
role of which is accomplished by llie weighting w^ 



(23) 



333 


r W, 


L<V V J 2w t 


G' » 


(24) 


12 




(25) 


where P is proportional to Opik's (64] earth 
encounter probability, and is normalized from the 
values which McCrosky and Posen [18] tabulated 
for Whipple's [42] cosmic weight w by 



The lunar weighting w^ from equation (24) is 

applicable at least for the back side of the moon. 

The extent to which it may be applicable for the 
front side is not certain. Several authors have 
reported that the decrease in flux due to the moon's 
lower gravitational field, i. e. , to about 50 percent 
by equation (10) , is somewhat compensated on the 
earthward side by an increase in flux due to the 
focusing effect of the earth. In 1959 Beard [65] 
gave a factor of 2 as die maximum increase in 
impact frequency on the earthward side. Shoemaker, 
Hackman, and Eggleton [66] reported in 1962 that the 
maximum local increase in the lunar influx may be 
greater than a factor 2 but must be less than a factor 
3. In 1964 Hale and Wright [67] reported consider- 
able theoretical analysis of the problem. More 
recently Hartmann [68] wrote that "... the net 
effect on the earthward side of the decrease in flux 
due to die moon's lower gravitational field and the 
increase due to the focusing effect of the earth is a 
decrease in flux by 0. 8, consistent with Opik [69] . . . " 
But Opik's [69] value 0. 80 is just the reciprocal 
of his value 1. 3 for the factor by which the influx of 
meteoroids into the atmosphere would be gravitation- 
ally enhanced if the geocentric velocity V^, were 

uniformly 20 km/sec. Convenient formulation for 
this calculation was shown by Shoemaker, Hackman, 
and Eggleton [66] . It may very well be that the 
lunar weighting w^ is just as applicable for the 

earthward side as for the back side. 

Further familiarity with the role of the three 
weighting functions w , w^, and w g is gained by 

considering their extreme values. By equations 


(23) through (25) which define the functions, each 
function has a unit mean value for either of the two 
sub-samples "A" or "B. " Consequently, those 
meteors which are common to both sub-samples 
have relative weighting values somewhat different 
in "A" than they have in "B. " In sub-sample "A": 

0. 4 £ w^_ ^ 2. 4, with 6 values w^ < 0. 6 and 8 values 

w > 1. 3; 0 is w :£ 2. 3, with 6 values w^ < 0. 2 and 

7 values w, > 1.7; and 0 ^ w ^ 7. 0, with 28 values 
1 s 

w < 0. 1 and 8 values w > 3. 5. Thus, as should 
s s 

be expected, the terrestrial weighting w departs 

from uniformity only mildly, while the lunar weight- 
ing Wj is rather more pronounced, and the spatial 

weighting w g is severe enough to reduce appreciably 

the effective size of the sub-samples. 

Statistical, Comparative, and 
Analytical Results 

The numerical results of statistical analysis 
with mass values from "Model-C" versus "Model- 
B" for Opik's [2, 28] compact stone versus stone 
dustball meteoroids are tabulated for comparison 
in Figure 14. The compared values are in very 
near agreement, except that the mass values by the 
dustball "Model-B" have a larger standard deviation. 
No basis is evident for any suspicion that the results 
from the further statistical analysis could be 
invalidated by using only the dustball "Model-B" 
mass values to represent Opik's [2, 28] results in 
comparison to those from "Model-A. " 

The quantitative results of the statistical analysis 
with sub-samples "A" and "B" are tabulated for 
ready reference and comparison in Figures 15 through 
17. For example, Figure 15 gives the weighted 
mean velocity of air-entering and lunar (primary) 
meteoroids as 16. 17 ± 0. 17 and 13. 66 km/sec, 
respectively, for sub-sample "A, " and 18. 83 ± 0.32 
and 18. 04 km/sec, respectively, for sub-sample 
"B. " In Figure 16, for example, air-entry velocity 
V has the correlations 0. 020 and 0. 191 with the 

CO 

"Model-A" mass m^ in sub-sample "A" and with 

the "Model-B" mass m in sub-sample "B, " 

respectively. The corresponding correlations with 
the lunar impact velocity, as approximated by the 
weighted geocentric velocity V in Figure 17, are 

(j 

0. 035 and 0. 228, respectively. 


13 



The various numbers tabulated in Figures 15 
through 17 are quite intricately interrelated, as a 
typical example of these intricacies, consider the 
meteor beginning height H . Hawkins and Soutlnvorth 

[53] showed how 11 for meteors photographed 

ID 

like those in the present sample, is related to 
velocity V . Meteors for which V is about 40 

* OC cc 

km/'sec have values of H about 100 km, the 

.D 

standardized reference height in the specification 

of absolute magnitudes, e. g. , M . The mean 

Pg 

values for H from Figure 15, S2. 99 i 0.23 and 

So. 14 ± 0.31 Ion with the terrestrial weighting w 

for sub-samples "A" and "B, ’’ respectively, are 
a reflection of the lower mean values for velocity 

V , already cited. The same explanation accoiuits 

foi the higher mean values of H witli the lunar 

ID 

weighting w^, S4. S4 and 87.78 km for sub-samples 

"A" and "B, 11 respectively: i. e. . the air-entry 
velocity V would lend to be distributed more toward 

OO 

higher values, except for the gravitational enhance- 
ment of the influx of slower meteoroids. Similarly 
for the even higher mean values of 11^ with the 

spatial weighting w i.e. , particles with higher 

closing velocity V are less favored for interaction 

with the earth, for more reasons tnan just gravita- 
tional attraction; otherwise, the air-entry velocity 

V would tend toward even higher values. Higher 

V is also the explanation of higher 11 „ for ''B" 

co p, 

versus "A :I with each of the four weighting functions 
(Fig. 15). In Figures 1G and 17, the correlations 
between and most of the other parameters are 

due primarily to the correlations between those same 
parameters and velocity V . 

Instead of the elongaLion A of the meteor radiant 
from the earth-way apex, its longitudinal projection 
(along the ecliptic) A is used as one of the 12 main 

A 

parameters in this study. It is related to A and the 

radiant celestial latitude ft by 

e J 

cos A = (cos A)/cos , (27) 

wherein it must be noted that \ and A are always in 
die same quadrant because cos (i never changes sign. 


Figures 1G and 17 show that the longitude 
interval A^ is correlated directly with perihelion 

q and inversely with the velocities V and V,.. But 

1 jo G 

q is inversely correlated with V and Y_. The 

CO O 

mean values of V and V from Figure 15 are lower 
co G 

for sub-sample ''A’ 1 than for "B, " and lower for 
the weighting functions w than for w^, which are 

lower than for w . These velocity changes are the 
s 

basis for some trends in the distribution of A in 

A 

Figures 18 through 20. Those figures show that, 

for either of the two sub-samples with either of the 

throe weighting functions, A has only one mode 

A 

(point of steepest ascent of the relative cumulative 
distribution) , always near the median (point where 
the relative cumulative distribution is 0. 5). In the 
respective Figures (18 through 20), by progressing 
from to Wj to w g in parts "A" or parts "B, 11 

or from part "A" to part ''B" ill either figure, the 
mode shifts progressively toward lower values, 
becoming very nearly 90 degrees for the sub-sample 
"B" with spatial weighting w . The combination 

"B" with w , giving A the mode near 90 degrees 
s A 

(Fig. 20b). is seen, in Figure 15. to represent 

a hypothetical sample with 26. 39 km/ sec average 

V . But the combination "A' 1 with w, . giving A 
co t A 

the mode near 105 degrees (Fig. 18a) , corresponds 
to 16. G2 km/ sec average V (Fig. 15). It could 

be expected that, by extrapolation to a sample with 
even higher values of V^, the mode of A^ might 

be shifted downward from 90, perhaps to as low as 

GO degrees. Davies and Gill (GO] found a A mode 

A. 

of about GO degrees with radio meteor data. 

The weighted means of meteor heliocentric orbit 
eccentricity e, tabulated in Figure 15, are within 
0. 02 of the corresponding median values, which are 
determined to about that accuracy by smoothing the 
cumulative distribution plots in Figures 21 through 
23. Because e is strongly correlated with velocity 
(Figs. 16 and 17), a progression of changes, from 
sub-sample to sub-sample and from weighting 
function to weighting function, is evident in Figures 
21 through 23. For example, with either "A" or "B, " 
the proportion of meteoroids in the total influx with 
e £= 0.25 is about twice as large for the earth (Fig. 

21) as for the moon (Fig. 22). 
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The medians of the weighted values of the 
perihelions a of the heliocentric orbits of the 
meteoroids are appreciably different from the 
corresponding mean values from Figure 15. With 
the terrestrial weighting w , the medians of q for 

sub-samples "A" and "B" in Figure 24 are 0. 909 
and 0. 8S1 AU, respectively, while the corresponding- 
means in Figure 15 are 0. 860G ± 0. 0054 and 0.8161 
t 0.0071 AU, respectively. With the lunar weighting 
w^, the medians of q for "A" and "B" in Figure 25 

are 0.886 and 0. 838 AU, respectively, while the 
corresponding means in Figure 15 are 0. 8427 and 
0.7804 AU, respectively. And, with the spatial 
weighting w , the medians of q for "A" and "B" 

in Figure 26 are 0. 844 and 0. 775 AU, respectively, 
while the corresponding means in Figure 15 are 
0. 8184 and 0. 7103 AU, respectively. The inverse 
correlation between perihelion q and the velocities 
V and V_ is evident both by visual inspection of 

CO (j 

Figures 24 through 26 and by comparison of the 

figures just mentioned for the means or the medians 

of q from "A" to "B" and from w, to w, to w . 

t 1 s 

The distributions of meteoroid perihelion q and 
aphelion q' with lunar weighting w^ and with spatial 
weighting w are of particular interest. They con- 
cern the variation of particle concentration with 
heliocentric distance and the corresponding impact 
flux against a moving interplanetary vehicle. In the 
data for "A" and "B, " q is given to only the nearest 
hundredth AU, and q' to only the nearest tenth. 

This grouping results in some irregularity of the 
cumulative distributions when several members 
have the same indicated values. Histograms of 
the probability density function corresponding to 
these 0. 01 increments as q approaches unity in 
Figures 25 and 26 are plotted in Figure 27. The 
histograms for the four largest abscissa values 
(0. 9C t) s l. 00 AU) are fitted in Figure 27 with 
dotted lines rising appreciably as q increases. This 
effect is more accentuated with "A" rather than "B" 

and with w, rather than w : i. e. , 

1 s 


f (q) = 145. 5 q - 137. 5 
= 64. 5 q - 59. 5 
= 90. 0 q - 83. 8 
= 21. 5 q - 18. 3 


for "A" and Wj \ 

for "A" and w / 

S 1(28) 

for "B" and w^ ( 

for "B" and w ) 

s ' 1 


Nevertheless, at lower values of q the cumulative 
distributions in Figures 25 and 26 are nearly linear 
over a wide range, indicating a constant value of 
f(q). Over the linear segments in Figures 25 and 
26, the relative cumulative distribution spans from 
0.25 to 0. 62 sloping 2. 00 for "A' 1 with w^, from 

0.36 to 0.70 sloping 1.86 for "B' : with w , from 

0.25 to 0.75 sloping 2. 38 for "A" with w , and from 

0.30 to 0.75 sloping 1.72 for "B" with w_. The 

nearly constant slope of these segments, which may 
measure f(q) more nearly free from selection 
effects, are extended by dashed horizontal lines 
into the plots in Figure 27. 


The distribution of meteoroid aphelion q' for 
different sub-samples "A" and "B" and different 
weighting functions w , w^, and w arc plotted in 

Figures 28 through 30. In each case the distribution 
is highly skewed toward the lower values of q', the 
most prominent mode is at 1 AU, and there is a 
slight secondary mode near 4 AU. With the spatial 
weighting w in Figure 30, there is another mode 

near 5 AU (Jupiter's heliocentric distance; . This 
mode near 5 AU in Figure 30 is barely evident for 
"A,” but for "B" it is accentuated about the same 
as the mode at i AU. Only values of q' less than 
10 AU are plotted, but the off-scale values and the 
"hyperbolics" (those known to have large but maybe 
finite q') were also used in the computation of the 
relative cumulative distribution. This is why the 
plotted values do not span the entire ordinate from 
0 to 1. Due to skewness, the mean values of q' 
from Figure 15 are all much larger than the 
corresponding medians from Figures 28 through 30. 
Nevertheless, the trend toward larger q' is quite 
evident for both means and percentiles in the 
aforementioned transition, "A" to "B" or w to w 
to w . t 1 

s 

A histogram of the probability density function 
of meteoroid aphelion q', corresponding to the 0. 1 
increments of q' near unity in Figures 29 and 30, 
is plotted in Figure 31. The same analysis is given 
for q' in Figure 31 as was given for q in Figure 27. 
The histograms for the four smallest abscissa 
values (1. 0 s q' < 1.4 AU) are fitted by linear 
segments; i. e. , 
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f(q') = - 0.48 q' + 1. 02 
= 0. 14 q' + 0.22 

= - 0.34 q' + 0.75 
= 0. 30 q' - 0. 00 


for "A" and w^ \ 

for "A" and w / 

S ) (29) 

for "B" and w ( 

for "B" and \v . I 
s 


For values of q beyond 2. 0 AU in Figures 29 and 30, 
there are nearly linear segments of the relative 
cumulative distribution, with "A" in the interval 
2 < q 1 < 4 AU, sloping 0. 19 and 0. 17 for w^ and w g , 

respectively, and with "B" in the interval 2 *' q' 

£ 5 AU, sloping 0. 15 and 0. 12 for w^ and w , 

respectively. It is not clear whether the slopes of 
these segments (which are extended with dashed 
horizontal lines into Figure 31) may represent 
f(q') as accurately as equation (29) in the interval 
l,0i q' 5 1.4 AU. It seems likely that equation (29) 
for f(q') should be a more substantial result than 
equation (28) for f(q) because the data spanned an 
increment of 0. 4 instead of 0. 04 AU. 


concentration. The velocity and volume vary as 
_ 1 

r 2 and r 3 , respectively, (e.g. , Opik [70] formu- 
lated the closing velocity V as the product of the 

_ 1 

earth's orbital velocity, proportional to r 2 , and a 
function of the meteoroid orbit parameters) . In 
interpreting the zodiacal light, it is usually assumed 
(e.g., by Fesenkov [ 7 1 J , Ingham [72], and Beard 
1 73] ) that the concentration of light-scattering dust 
particles varies with r _n , where n is a constant. 

In this study a similar assumption is made for the 
flux F of particles in the vicinity of an interplanetary 
vehicle at heliocentric distance r. Then, the value 


of -n in the corresponding exponential r n l is 

1 1 
found by subtracting 3 for volume and — for velocity 

from equation (30) ; i. e. , 


f ~ r £ W - f (q') " 3 - 5 = r "l 


nr i 

(l/F) — — = -n, r _1 = -n, 1 AU 

dr 1 1 


Rejecting equation (28) and using the dashed 
lines for f(q) in Figure 27 while using equation (29) 
for f(q') gives the expected values for differences 
between the probability density functions of perihelion 
q and aphelion q' when both q and q' are in the 
vicinity of 1 AU; i. e. , 


- n^ = f(q) - f(q') - 3. 5 with w^ I 

= -2.04 for "A" / (32) 

= -2. 05 for "B" . | 


f (q) - f(q') = 1.40 
= 1.45 
= 2. 03 
= 1.48 


for "A" with w 
for "B" with w^ 
for "A" with w 

s 

for "B" with w 

s 



Equation (30) represents the derivative at 1 AU, 
with respect to heliocentric distance r in AU, of 
the relative number of meteoroids available with 
orbits appropriate for the constitution of a flux of 
particles in the vicinity of a vehicle which has an 
approximately circular heliocentric orbit. This is 
the reason for the lunar weighting w^, the moon is 

such a moving target with a practically circular 
heliocentric orbit which remains always approximately 
at 1 AU. The corresponding impact flux varies as 
the product of the closing velocity V and the volume 


A corresponding exponential ~ r s is assumed for 
the concentration of meteoroids in space, except 
that there is no component for velocity variation 
and equation (31) with w is used instead of equation 

(30) with w, i.e. , 


-n g = f(q) - f(q'i - 3 with w I 

= -0.97 for "A" > (33) 

= - 1. 52 for "B" 


It may not be particularly significant that, in 
comparison to zodiacal light and solar F-corona 
results, "A" in equation (33) agrees with Fesenkov's 
[71] result that particle concentration varies with 
r -1 , while "B" agrees with Ingham [72], and Beard 
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[73], that concentration varies with r- 3 / 2 . The 
linear dimensions of the meteoroids in "A" and "B" 
are perhaps two orders of magnitude larger than 
those most representative of the optical effects. 

Still, the dissipative optical effects are reported 
(e.g. , by Briggs [74]) to be quite pertinent to the 
distribution of meteor particles in space. 

In 1954 Whipple [42] noted that the orbits of 
comets can be distinguished from those of asteroids 
by the value of a parameter K, 

K = log (f^r) - 1 ’ 

where, as Jacchia and Whipple [75] noted sub- 
sequently, the term in parentheses is inversely 
proportional to the square of the velocity at aphelion. 
In suggesting the criterion, Whipple noted that K is 
positive for 96 percent of the known comet orbits 
and negative for all but 3 asteroid orbits. Jacchia 
and Whipple [75] cautioned that the values of K are 
not well known for asteroids of perihelion q < 1 AU. 
The practice by some authors (e. g. , Whipple [42] , 
Hawkins and Southworth [76], Jacchia and Whipple 
[75], Whipple and Hawkins [12], and McCrosky 
and Posen [ 18] ) of applying the criterion to the 
orbits of meteor particles has been criticized 
(e.g. , by Dycus and Bradford [77] and by Ceplecha 
[44]) as an inappropriate extrapolation (with respect 
to perihelion distance) from the data on which it was 
based. Still it is an interesting parameter; and the 
present analysis gives results quite different from 
others. Whipple [42] analyzed 144 meteors "from 
double-station photographs by cameras equipped 
with rotating shutters . . . presenting a fairly 
representative sampling of meteors to about 0 to 1 
magnitude visually. " His results, illustrated 
graphically by Whipple and Hawkins [12] and 
further discussed by McCrosky and Posen [18] , 
showed only 5. 6 percent with K < 0. In their own 
analysis of the 2529 photographic meteors, from 
which sub-samples "A" and "B" were taken for the 
present analysis, McCrosky and Posen [18] found 
25. 5 percent with K < 0. But in the present analysis 
with terrestrial weighting w^, the meteors with 

K < 0 in sub-samples "A" and "B" are 77 percent 
and 67 percent, respectively. That (essentially) 
meteors with lower than a limiting mass were not 
included in "A" and "B" does not explain the 
preponderance of members with K < 0, because 
Figure 16 shows significant (albeit small) positive 
correlation with respect to mass. With probable 


error less than 0. 037, the correlations between K 

and in. in "A" and between K and m„ in "B" are 
A B 

0. 084 and 0. 160, respectively, with terrestrial 

weighting w^. If one would interpret K < 0 as some 

indication for asteroidal origin of meteor particles, 

then the positive correlation with mass, especially 

in "B, " would support Kresak's [20] finding that 

the ratio of asteroidal to cometary meteors increases 

with decreasing brightness. And, except for Opik's 

[78] cautioning that not only dustballs but also 

meteoroids of compact iron and of compact stone 

may originate from comet nuclei, the positive 

correlation between the K parameter and mass 

also supports Ceplecha's [44] finding "that weaker 

meteors photographed by Super-Schmidt cameras 

contain more observed dense particles than the 

small camera meteors" (i.e. , the brighter meteors). 

But a further confirmation that the smaller meteoroids 

actually do have more structural integrity is given 

through the appreciable positive correlation between 

meteor beginning height H and mass. From Figure 
B 

16, the correlations between H„ and mass m. in 

B A 

"A" and mass m in "B" are 0. 131 and 0.249, 

B 

respectively (with weighting w ). This tends to 

imply that the smaller meteoroids are more likely 

to have the higher structural integrity to remain 

optically quiescent by delaying breakup down to the 

lower meteor beginning heights II . Jacchia, 

B 

Verniani, and Briggs [79] found, by analysis of 
413 precisely reduced Super-Schmidt photographic 
meteors, that on the average, meteors in short- 
period orbits appear to be about 1.4 times as dense 
as meteors in long-period orbits, and also, for 
short-period meteors, that the density increases as 
the aphelion q' decreases. Their results are 
supplemented by those of other authors (e. g. , 
Babadzhanov and Kramer [80]) who find that the 
proportion of short-period orbits is higher among 
the faint meteors than can be explained by observa- 
tional selection. Thus, the smaller orbits tend to 
have the smaller meteoroids, tending also to have 
the higher densities and (whether coincidentally or 
significantly) the lower (most asteroid-like) values 
of the criterion K. The cumulative distributions of 
K with the weightings w^, w^, and w are shown in 

Figures 32 through 34, respectively. In each case 

not only the correlation with mass but also the 

median of K is higher (more comet-like) for "B" 

than for "A," for w, than for w., and for w than 
~ its 

for w . 
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The density of meteor radiants per square 

degree of the celestial hemisphere is proportional 

to the slope of the relative cumulative distribution 

of sin 8 , as has been mentioned. Using \(1 j 
e e 

instead of some attention was given to the 

distribution of this parameter in the development of 
the weighting functions. The weighting function w 

gives a distribution of sin |/3 j which approximates 

that of sin p when only positive values of j3 ^ are 

counted and they are weighted as a function of the 
apparent fraction of the circle of celestial latitude 
through the meteor radiant. The results, plotted in 
Figure 10, show that half of the meteor radiants in 
the atmosphere are within ± sin -1 0.44 or 26° and 
± sin -1 0. 39 or 23° of the ecliptic for sub-samples 
"A" and "B, " respectively. Also, from the slope 
of the plot in Figure 10, one sees that between ± sin -1 
0. 32 or 19° and i sin -1 0. 90 or 64° the density of 
radiants per square degree is much less than else- 
where. It is of some interest to check these results 
with those with the lunar weighting w^. Figure 35 

shows that the intermediate values of S/3 ' are still 

e 

associated with the low values of radiant density per 
square degree. The median values of 1/3^ | are 

shifted to somewhat higher values with the lunar 
weighting w being sin -1 0.47 or 2S° and sin -1 0.45 

or 27° for "A" and "B, " respectively. The density 
of meteor radiants per square degree, as they enter 
the atmosphere or as they approach the moon or an 
interplanetary vehicle, is of interest for impact 
hazard considerations. Figures 10 and 3G show that 
the prospect for shielding a vehicle against an 
enhanced meteoroid flux from regions in the vicinity 
of the ecliptic are of little practical merit because 
even within angular distance from the celestial poles 
up to 26 degrees the radiant density per square 
degree is just as high as it is within ± 19° of the 
ecliptic. Anticipation of such practical shielding 
possibilities undoubtedly results from information 
about the distribution of sin i, where i is the 
inclination of the meteoroid orbit to the ecliptic. 

This distribution is plotted in Figure 36, where the 
spatial weighting w is used to approximate the 

characteristics of the meteoroids in space. The 

median values of sin i in Figure 36 are 0.24 or 

sin 14° and 0.29 or sin 17° for "A" and "B, " 

respectively. The mean values of sin i in Figure 

15 show a pronounced and consistent tendency to 

increase from"A" to "13" and from w, to w, to w . 

t 1 s 


The cumulative distributions of meteoroid air- 
entry velocity V (with terrestrial weighting w^) 

and lunar impact velocity V^, (with lunar weighting 

w.) are plotted in Figures 37 and 3S, respectively. 

The medians for sub-samples "A" and "B, ” 
practically coincident with the single distinct modes, 
are 15.3 and 16.5 km/sec, respectively, for V 

CO 

in Figure 37, and 12.6 and 14.8 km/sec, respectively, 
for V in Figure 3S. These medians are consis- 

tently lower than the corresponding means, tabulated 
in Figure 15, because the distributions are skewed 
toward the low velocity ends of the ranges of values. 

By multiplying the abscissa by r _1 / 2 j Figure 38 
would approximate the cumulative distribution of 
meteoroid closing velocity V with respect to an 

interplanetary vehicle in a (practically) circular 
orbit at heliocentric distance r Ali, 

V = V_ r" 1/2 . (34) 

V (j 

Because the vehicle orbital velocity varies with 
r" 1 / 2 , equation (34) follows from Opik's [70] formu- 
lation for closing velocity proportional to the product 
of orbital velocity and a function of the meteoroid 
orbit parameters, together with the assumption that 
the particular function of meteoroid orbit parameters 
does not depend strongly on r. 

In addition to the correlations between the air- 
entry velocity and the mass evaluations ni^ 

(in sub-sample "A") and m (in "B"), further 

13 

information about the relations between these 
parameters and the differing constitutions of "A" 
and "B" is plotted in Figure 39. Each of the two 
sub-samples "A" and "B" was ordered with respect 
to its characteristic mass parameter (m for "A" 

and m for "B") . Each of "A" and "B" was 

J3 

partitioned (by au intermediate value of its mass 
parameter) into larger-mass and smaller-mass 
sets giving the most nearly balancing sums of the 
terrestrial weighting w,. The relative cumulative 

L 

distributions of V for the different mass regimes 

OO 

are shown in Figure 39 to be significantly different 

for "B, " but not for "A." This checks with the 

corresponding correlation values from Figure 16; 

i. e. , 0.020 with m in "A" (less than the probable 
A 


IS 



error of 0.037) and 0.191 with m_ in "B. " The 

a 

meteors in Figure 39 with small values of the mass 
(dot points) tend to have the lower values of 

velocity in "B" than the meteors with large 

values of m^ (asterisk points) . 

In making statistical studies with meteor data, 
most investigators (e.g. , Hawkins and Upton [81], 
Spik |32], Lindblad [82], Watson [83], Lovell [4, 
84), McCrosky [85], Whipple [86], Dohnanyi [87], 
and Erickson [88]) use the distribution of absolute 
magnitude in a random sample of meteors for the 
determination of the mass distribution. Their 
results are based on the experimental finding that 
the "luminosity function" N > ^ (M) , the number of 

meteors equal to or brighter than a given absolute 
magnitude M (photographic M or visual M ), 
is of the form 


n >m (M) = n >m (0) * 


M 


(35) 


in a random sample of meteors, where the parameter 
X (photographic Xp^ °r visual x y ) is considered to 


be constant over a fairly wide range of magnitudes 
M. Then, by representing magnitude M as a 
linear function of log mass m (e.g., m or m ) 

A jj 


at constant velocity, etc. , 


M ( m) = Iq, + log + k 2 log m 


+ k 3 log cos Z R 


(36) 


where kj, k 2 , and k 3 are constants, it follows 

(as Lindblad [82] shows) that the "mass function" 

N (m) , the number of meteors of mass equal to 
>m 

or greater than m, is 


investigators who have followed that procedure have 
taken the theoretical value -5/2 for k 2 (i.e. , "5" 
and "2" because when magnitude is increased by an 
increment "5" the luminous intensity is increased 


by a factor 10 


) . In an evaluation of x in 
Pg 


equation (35) for McCrosky and Posen's [18] entire 
sample of 2530 meteors, Dohnanyi [87] used the 
empirical value -2. 25 for k 2 in equation (36) which 
Jacchia, Verniani, and Briggs [56] found by statisti- 
cal analysis of photographic meteor data. Also, 
corrections which Dohnyani [87] made for the 
photographic magnitude scale discrepancy, noted 
by Kresak [20] for the 300 meteors tabulated by 
Hawkins and Soutliworth [53] and also, with some- 
what different magnitudes M by McCrosky and 

Posen [18] , have the effect of reducing the value 
of -k 2 by a factor 0. 79 for use with McCrosky and 
Posen's [IS] data. Then the effective value of k 2 
in equation (36) should be -1. 78 in the luminosity 
function for present analysis with sub-samples "A" 
and "B. " 


The present analysis, with data selected by a 

criterion for absolute magnitude M extrapolated 

Pg 

for a definite limit of velocity V , is effectively 


more nearly selected down to a definite limit of 
mass. Consequently, the distribution of mass is 
determined directly, instead of from the luminosity 
function. The distribution of meteoroid mass in the 
terrestrial and lunar influxes and in the concentration 
of particles in interplanetary space are plotted in 
Figures 40 through 42, respectively. The straight 
lines approximating the respective plots were not 
fitted by least squares, but were determined by the 
points for which the ordinate, the log relative 
cumulative distribution, has the values -0.25 and 
-1. 00. These two points are appropriate for sub- 
samples with 333 members because the upper point 
is near the region beyond which the log cumulative 
distribution has what must be a spurious curvature, 
and the lower point is near the region where the 
slope is unduly influenced by the randomness of the 
individual data points. 


N (m) = N (1) m 
>m >m 


k 2 log x 


= N (1) m 
>m 


-(s - 1) 


I 

\ 


(37) 


where s, following the convention of some authors 
(e. g. , Divari [89]), is a constant over a fairly wide 
range of in. Except Dohnanyi [87], all of the 


That the indicated lines in Figures 40 through 
42 must be the best models for the populations of 
meteoroids, which are represented by the sub- 
samples "A" and "B" statistically weighted with 
w , w R or w g , is not evident beyond a reasonable 

doubt. But no process or procedure is evident by 
which the results could be interpreted better when 
any one of the plots is considered, whether individ- 
ually or in comparison with another. The slopes of 
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the plots in Figures 40 through 42 (involving the 
respective weighting functions w^, w^, and w g ) , 

are -1.48, -1.42, and -1.42, respectively, with 
sub-sample "A" and are -1. 31, -1.21, and -1.37, 
respectively, with "B. " The slopes from Figure 
40, -1.48 and -1.31 for "A" and "B, " respectively, 
can be equated to the exponent k 2 log >; (T in equation 


(37) to find the values of the parameter x 


Pg 


• By 


using the empirical value 2.25 for k 2 in order to 
compare results on the same basis with most 
reported results, one finds 


X = 4. 55 inferred from "A" 
Pg 

= 3. 82 inferred from "B" . 


It must be understood that the values of x in 

Pg 

equation (38) are inferred for a random sample of 
meteors, like McCrosky and Posen's [18] entire 
sample of sporadic meteors instead of the selected 
sub-samples "A" and "B, " and that they presuppose 
the statistical independence of the mass and velocity 
of meteor particles. Due to the correlation between 
mass and velocity (tabulated in Figure 16) it should 

be expected that the values of x determined 

Pg 

directly from photographic meteor data may be some- 
what different from those in equation (38) . Never- 
theless, the value 3. 82 for "B" in equation (38) 
should be considered a fairly close agreement with 
the value 3.44 which Hawkins and Upton [81] found 
directly for x p ^ in an analysis of the distribution 

of magnitude for 2123 sporadic photographic meteors. 


The agreement would be even almost perfect if the 
theoretical value -5/2 were used for k 2 . 

In 1941, and in a revised edition in 1956, Watson 

[S3] gave visual meteor data showing that in the 

interval -3 £ M £ 2 the number of meteors seems 
v 

to increase approximately 2. 5 times with each 
fainter magnitude. Watson [83] interpreted that 
result to indicate that log x ' is exactly the reciprocal 

of -k 2 in equation (37) , so that the value of the 
corresponding parameter s is 2. Many authors 
(e.g. , Divari [89]) have quoted Watson's [83] 
results in that sense. But the same value of k 2 
cannot be applied to both visual and photographic 
meteor data. The visual data are complicated by 
the Purkinje effect (transition from phototopic to 
scotopic visual sensitivity, i. e. , from cone to rod 
sensors) and other differences between visual and 
photographic spectral sensitivity which are corrected 
collectively by a color index C, 

C = M - M . (39) 

pg v 


Verniani [26] suggests that, in the interval 

-2 < M £1.5, -Jacchia's [90] results for color 
Pg 

index C can be approximately described by 


C = 0. 28 M -1.34. 
Pg 


(40) 


By equations (39) and (40) the derivative 

with respect to M is 0.72. Therefore, 
/am 2 Pg 


of M 

v 

by equation 


2. For meteors sufficiently bright to be effectively free of both visual and photographic selection effects, 

equation (35) gives the same number of meteors with subscripts for visual or with subscripts for photographic 

meteors. Then the relative derivatives of these numbers N and N with respect to M are also equal: 

, _ v pg pg 



dN 

ES. 

dM 

Pg 


log x 
e 'pg 



dM 

(l0g e V 


Pg 


0.72 log e X v 


l0g e 


X 


0.72 

V 


X 


Pg 


= X 


0.72 

V 
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X 


pg 


0.72 for _ 0 _ 1 < M < 2. 4 . 
V V 


(41) 


Watson's [ 83 ] value 0.4 for log in equation 

(41) would give 0. 288 for log x ^5 and, with 

Pg 

Jaceia, Verniani, and Briggs' (56] value -2. 25 for 
k 2 from equation (36) , the corresponding value of 
the mass distribution parameter s in equation (37) 
would be 1. 65 instead of 2 over part of tire range 

M in Watson's [83] data. This reduces consid- 
v 

erably the basis for Watson's [83] inference about 
the mass distribution. The somewhat higher value 
3. 7 for x v was reported in 1957 by Millman and 

Burland (see Millman [91]) for a larger sample of 
visual sporadic and shower meteors. In 1965 
Lindblad [82] reported an observational program and 
analysis of visual meteors directed toward the ident- 
ification and elimination of tire major sources of 
bias in visual meteor data. From a total of 4000 
meteors, 859 were determined to be sporadic and 
255 of those sporadic were brighter than the faintest 
absolute magnitude at which the survey was judged 
to be complete; i. e. , free of selection effects. The 
value of the parameter x y was reported to depend 

on the arbitrary width of that part of the observers' 
field of view within which an observed meteor would 
bo counted. Extrapolation of the sporadic meteor 
data to zero angular distance from the center of the 
field of view gave 4. 3 for the parameter x y - The 

difference between Lindblad's [82] value 4.3 for 
X v and Hawkins and Upton's [81] value 3.44 for 



role of the color index over part of the involved 
range of visual magnitude. 


It was noted above that the value 3. 82 for x , 

Pg 

which by sub-sample "B" one would infer for a random 
sample of meteors, agrees well with the value 3.44 
found directly by Hawkins and Upton [81] , and that 
the theoretical value -5/2 for k 2 would give 3.34 

for x • This does not seem to be a sufficient basis 
Pg 

to question the merit of Jacchia, Verniani, and 
Briggs’ [56] statistical determination -2. 25 for k 2 . 

On the contrary, if one takes, as Erickson [88] did, 
the (above) experimental value -2. 25 for k 2 , then 
the Hawkins and Upton [81] value 3.44 for x 

Pg 

would imply a value of 2.21 instead of 2.34 for the 
mass distribution parameter s in equation (37) if 
mass and velocity could be assumed independent. 


Actually, the value 2. 21 is not s but the correspond- 
ing parameter s g in the distribution of kinetic 

energy for meteor particles; i. e. , 


N 


> A mV 2 


UmV^ 


• imV‘ 


(l)(4mV^) 


-(s 


- 1 ) 
(42) 


The distribution of meteor particle kinetic energy, 
unlike either mass or momentum, is practically 
well represented by a random sample of meteors. 


Fialko [92] showed that the distribution param- 
eters s for mass m in equation (37) and s^ for 

energy' in equation (42) would be identical except 
for any statistical dependence between the air-entry 
mass and velocity. The English translation of 
Fialko's [92] paper from Russian appeared in the 
September-October 1965 issue of Soviet Astronomy- 
AJ. Being unaware of Fialko's [92] analysis, the 
author had distributed a Trip Report in August 1965 
with a short proof of the more general result that, 
in equation (37) for the distribution of particle mass, 
the parameter s is invariant with respect to the 
replacement of mass by the product of mass and any 
exponential function of velocity, provided only that 
mass and velocity were statistically independent. The 
same proof was also published by Dalton [93] in a 
NASA Technical Memorandum in November 1965. 
Marcus' [94] boundary-effect criticism of Fialko's 
[92] analysis would apply equally to Dalton's [93] 
more general result, but Fialko [95] showed that 
those boundary effects are not practically significant 
for meteoroids in the terrestrial influx. Subsequently, 
in an analysis of 184 radio meteors, Fialko [96] 
found a value & 2. 3 for the energy distribution param- 
eter s^ inequation (42). But, as Fialko [96] noted, 

"quite apart from the complexity of this method, one 
should realize that so far an accurate dependence of 
the meteor ionization 'probability' p on V has not 


been obtained. 1 


The value 2.21 for s inferred from 
e 


the Hawkins and Upton [81] value 3.44 for x is 

Pg 

probably the most substantial result for the energy 
distribution, at least for sub-sample "A. " 


With sub-sample "B" the correlation 0. 191 
between mass m and velocity with terrestrial 

weighting w might be sufficient to account for the 

difference between the value 2.31 for tire parameter 
s in the mass distribution and the value 2. 21 for the 
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corresponding parameter in the energy distribu- 
tion. But the corresponding results for "A, " with 
the value 2. 48 for s and only 0. 020 for the mass- 
velocity correlation, are somewhat less plausible. 

But in either case "A" or "B, " as long as one has not 
yet made a choice between "A" and "B, " the best 
basis at hand for the computation of the meteoroid 
impact hazard in space is the value 2.21 for the 
parameter in the energy distribution, equation 

(42) , for both the terrestrial influx and the lunar 
influx. Near the end of the last part of this section 
a choice is made between "Model-A" and "Model-B. '' 
There some further discussion of the role of k 2 and 
its effect on s^ is g'iven. In constructing the models 

for the distributions of energy, it is noted that the 

lines representing the distributions of mass in 

Figures 40 and 41 were fitted through the points for 

which the ordinates were -0.25 and -1. 00. The 

region of best fit tends to be centered approximately 

around the quartile, where the ordinate is -0. 60. 

Through the point with this ordinate -0. 60 on the fitted 

line, with abscissa m , a line with slope -1.21 for 
e 

-(s - 1) has the ordinate intercept 


the lunar weighting w^ the expected values of V^ 42 

are (14. 72) 2-42 and (22. 08) 2 - 42 for "A" and "B, " 
respectively. 

For "A" and "B" in Figures 40 and 41, the mean 
of the four values of s for mass is 2. 35. The effect 
of the correlation between mass and velocity is that 
the value of for energy is reduced to 2. 21. By 

linear interpolation (with respect to n in mV 11 ) the 
corresponding parameter s for momentum is 

2.28, and the ordinates y and y corresponding 
m in 

m 

to y and y , respectively, in equations (43) and 
e 

(45) , respectively, are 


y = - 0. 60 + 1. 28 log in 
in e 

m 

"333 

l V 1 ' 28 w 

v = y + log — 

•All -All H 333 

m v 


(46) 


(47) 


y = -0. 60 + 1.21 log m . (43) 

m e 

e 


This is the line for the mass distribution which would 
be inferred from the energy distribution if one could 
assume that mass and velocity were statistically 
independent. By representing the corresponding 
ordinate for the (similarly plotted) energy distribu- 
tion by y , Fialko's [92] formulation gives 


y = y + log E 
4 e -An ° 
e 


(V 2 /2)* e 


(44) 


where E[ ] designates expected value; i. e. , 


y 


e 



0. 36 + log- 


333 

Y V ’-« 

1 


333 



w 


(45) 


With the terrestrial weighting w the expected values 
of V^ 28 in equation (47) are (16. 69) 1-28 and (20. 71) 1,28 

for sub-samples "A" and "B, " respectively; and with 

the lunar weighting w the expected valuesof V 1 - 28 
1 G 

are ( 12. 55) 1-28 and ( 17. 38) 4 - 28 , respectively. 

In the same paper in which he introduced his 
comet-asteroid criterion K in 1954, Whipple [42] 
suggested that a weighting factor, V* 2 for V > 19 

OO OO 

km/sec and 3 V -2 for V < 19 km/ sec should reduce 

CO OO 

the distribution of observed velocities ill a random 
sample of meteors roughly to a distribution applying 
to meteoroids of a constant mass. He defined his 
"cosmic weight" w as 


w ~ P' 1 V -2 

w =° 


(48) 


With the terrestrial weighting w the expected values 
of V^ 42 in equation (45) are (17. 65) 2 - 42 and (23. 91) 2 ' 42 
for sub-samples "A" and "B, " respectively; and with 


where P is the meteoroid's earth-encounter probabil- 
ity which was formulated by Opik [64]. This weight- 
ing w is often tabulated (e.g. , by Hawkins and 

Southworth [76], Jacchia and Whipple [75], and 
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McCrosky and Posen [18]) and is of some interest 

to check with the present analysis. The weighting 

w was formed as the product of several factors, 
w 

one being V 2 because the effective surveillance 

area is considered to vary with the square of the 
meteor beginning height H , which in turn was 

Jo 

1 

considered to vary with V 4 . Thus, when meteor 

CO 

height is accounted for by a separate weighting, and 
one is concerned with the terrestrial influx instead 
of the total population of particles in interplanetary 
space, the weighting with respect to velocity would 
be ~ V^ 3 / 2 . Recently, Miller [51] derived a 

corresponding exponential velocity weighting for 
random photographic meteors. He found that the 
most appropriate velocity exponent is -2. 87 (s - 1) , 
where s is the parameter in the mass distribution, 
equation (37). The values of s corresponding to 
Figure 40 are 2.48 and 2. 31 for sub-samples "A" 
and "B, " respectively, and the corresponding values 
of the velocity exponent for Miller's [57] weighting 
are -4. 25 and -3. 76, respectively. An independent 
check on the exponent for velocity weighting of ran- 
dom meteors is afforded by assuming different values 
for die exponent and calculating the weighted mean of 
with a random sample. Then by interpolation 

one can find the value of the exponent which gives die 
same mean value of as is found in the present 

analysis with the selected sub-samples "A" and "B. " 
Tliis was done by hand computation with Hawkins and 
Soudiwordi's [53] histogram of velocity values 
classified in 5 km/sec intervals for a random 
sample of 285 sporadic photographic meteors. With 
uniform weighting except for velocity, and with 
velocity weighting exponents -2, -3, and -4 the 
mean values of velocity V are 20. 42, 17.48, and 

15. 77 km/sec, respectively. From Figure 15 the 
means of with uniform weighting are 16. 62 and 

20.23 km/sec for "A" and "B, " respectively. By 
interpolation one finds that the equivalent velocity 
weightings, without any adjustment for effective 
surveillance area, are ~ V -3 ' 50 and ~ V -2 ' 06 for 

0O OO 

sub-samples "A" and "B, " respectively. 


Metric of Cumulative Flux for 
Particle Parameters 


One can calculate the mean flux (or the expected 
future flux for the year-around average exposure) 
instead of just the relative rate. The 333 meteors 
in each of the sub-samples "A" and "B" were, in 
each case, all of die sporadic meteors in McCrosky 
and Posen's [18] tabulation which were bright 
enough for selection with respect to a criterion 
function of magnitude and velocity roughly correspond- 
ing to meteoroid mass. McCrosky and Posen [18] 
noted tiiat the 2529 meteors they tabulated were 70. 4 
percent of all meteors photographed during 1125 
hours of exposure, that 13. 8 percent had been 
rejected because the correspondence of points in the 
photographs from the two separated stations was 
uncertain, that 0. 6 percent were rejected due to an 
ambiguity of radiant, that 4. 8 percent could not be 
reduced for random technical reasons, and that 10.4 
percent were just too faint. The 10. 4 percent of 
too faint meteors are of no consequence to the 
present analysis because they would not have passed 
the selection criterion; but the other 19. 2 percent 
must be accounted for by reducing the 1125 hours of 
exposure time. Of the 2529 meteors which they 
tabulated, they designated 2066 sporadic and 463 
stream meteors. But only 2040 of the sporadics had 

values tabulated for M and m. , i. e. , could be 
Pg A 

considered for selection. Then, the effective total 
exposure t was 


t = 


1 - 


19,2 

70.4+ 19.2 


/ 2040 \ 
\ 2066 ) 


1125 hrs = 873 hrs. 


The effective surveillance area of equipment was 

given by Hawkins and Upton [81] as 5980 km 2 at 

the 90-km height for which the cameras are adjusted 

to overlap. The means H of the beginning heights 
13 

H (from Figure 15) for sub-sample "A" are 82. 99 
13 

and 84. 84 km for terrestrial weighting w and lunar 
weighting w respectively, and for "B" are 85. 14 
and 87.78 km for w and Wj, respectively. The 
effective area A is considered to be (H^/90) 2 2980 
km 2 . For both "A" and "B" the effective number N 

e 

of sporadic meteors, with respect to which the 
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indicated relative cumulative distributions must 
relate, is 333 with the terrestrial weighting w^; but 

with the lunar weighting' w the number of sporadic 

meteors which relate to the indicated area A and 
time interval t must be reduced by the factor 

from equation (10) , 0.47 for "A" and 0. 54 for "B. " 
But, in showing the cumulative flux with respect to 
mass, momentum, and energy, in specific units, 
the interests of most readers is probably served best 
by increasing the sporadic flux by the factor 2529/ 
2066 by which the total sample exceeds the sporadic 
component in the data shown by McCrosky and Posen 
[18]. Therefore, the increment Ay which must be 
added to the log relative cumulative distribution of 
sporadic meteors to get the log mean total flux per 
square meter per second is 


\ / 2b>9 \ 



B V 


)U>66 ) 

(333 k f ) 

10 c (19S0)I - 

-) (3GO0H873J 

( 


= - 9.453 ♦ log k* - .Mo g 7 i b . 

where 6 has the values 0 and 1 with w and w^, 

respectively. The results for Ay in equation (49) 
are -13. 30 and -13. 62 with weighting w and w^, 

respectively. The results for "A" and "B" are 

averaged and do not differ from their mean result 

by more than ±0. 015 with either w or w, . 

J t 1 

The mass cumulative influx F , found by 

>m 

adding Ay from equation (49) to the equations 
indicated on the plots in Figures 40 and 41, is 


log F = - 1. 48 log m . - 15. 05 
° >m A 

for "A" with w 
= - 1.42 log m^ - 15.29 
for "A" with w^ 

= - 1.31 log m B - 14.77 
for "B" with w 



= - 1.21 log m - 14. 93 
13 

for "B" with Wj 


(50) 

( con I’d) 


Also, in each case, the equation indicated on the plot 
in Figures 40 and 41 can be solved for the value 
log m of the abscissa corresponding to the ordinate 

value -0. 60. By substituting the results for log m^ 

into equation (46) one gets, for the mass reference 

intercept y for momentum, the values -1. 59 
J m 

m 

and -1.56 with w and w^, respectively, for "A," 

and the values -1.45 and -1.35 with w and w^, 

respectively, for "B. " Then, by equations (47) 
and (49) , the momentum cumulative influx F 

' ' ’ simr 


log F = -1. 2S log m. V - 13. 33 
>mv A °o 

for "A" with w 

= - 1.28 log m V - 13.78 
A G 

for "A" with w 

= - 1.28 log niV -13.07 

for "B" with w 

= - 1. 28 log mV - 13. 3S 
B G 

for "B" witli Wj 


(51) 


Similarly with the mass reference intercept y 

e 

in equation (43) , one gets the energy intercept y 

in equation (45), and adds Ay from equation (49) 

to get the energy cumulative influx F , i. e. , 

>imv 
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l 0 gF >imV 2 — 1 - 21 l 0 g( ^ 1 A V ^ 
for "A" with w 

= - 1.21 log ( 2 ln A V q) 

for "A" with w^ 

= - 1.21 log (im^) 

for "B" with w 

- - 1.21 log (i,n B V 2 G ) 

for "B" with w 



It is of some interest to compare the mass 
cumulative terrestrial influx in equation (50) with 
the results of tlie meteoroid-puncture experiments 
in near-earth-orbit space from Explorers XVI and 
XXIII and the three Pegasus flights. Dalton [16] 
showed that the product m 19 / 54 V 2 / 3 must be constant 
for just-puncturing meteoroids when the other impact 
circumstances are fixed. By taking 26.7 km/ sec 
for the mean impact velocity, Dalton [16] showed 
that (lie combined puncture data for the Explorers 
XVI and XXIII would indicate 10 _s - 374 impacts per 
square meter per second of particles of mass down 
to 10 -8 ' 17 gram, and that the corresponding values 
implicated by the Pegasus' thickest target would be 
J0-7.1G4 and to -6-03 , respectively, in the same units. 
With the mean values 16. 17 and 18. 83 km/ sec for 
"A” and "B, " respectively, (for V weighted w 

from Figure 15) the indicated just-puncture mass 
values should be increased by factors 10°' 413 and 
10 0-287 for "A" and "B, " respectively. Actually, 
this explanation is overly simplified because the 
E[V 2-42 ] rather than £[V ] is pertinent here, as 
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in equations (44) and (45) , but these values are not 
known for the minuscule particles. The results are 
plotted in Figure 43 with an extrapolation of the 
equation (50) mass cumulative terrestrial influx 
for comparison. But the mass values for the just- 
puncturing meteoroids plotted in Figure 43 pre- 
suppose the velocity distributions from the meteor 
sub-samples "A" and "B" to be applicable for 
particles of nearly eigiit orders of magnitude smaller 
mass. This presupposition is more nearly a valid 
representation of the results from "A" than from 


"B. " With "A" the 0. 020 correlation between mass 
m^ and velocity V^ is negligible in comparison 

with tlie probable error of ± 0. 037, but with "B" 

the 0. 191 correlation between mass m and 

JB 

velocity V is unmistakably significant. This 

would imply from "A" a lower velocity, and con- 
sequently higher mass for the satellite-puncturing 
particles. One can show that the Explorer puncture 
data most probably concerns typical zodiacal light 
particles. 

According to Takaknbo (see Ericke [57]) the 
radius of the dust particles most effective in the 
zodiacal light and the solar F-corona is about 20 
microns. With Opik's [2] value 3.4 g cm -3 for the 
average density of meteoritic stone, the mass of 
such a particle is 1. 14 x 10 -7 g. Other authors 
(e.g. , Briggs [74]) adopt density values as low as 
0.2 g cm -3 for zodiacal light particles of 20 micron 
radius. So the mass probably is between 10 -8 
and 1 0 -7 g. According to Opik [97] the velocities 
of the zodiacal light particles relative to the moving 
earth, the moon, or an interplanetary vehicle 
(i.e. , Vq) are small, the average being about 

one-fifth of the Earth's orbital velocity, or 6 km/sec. 

If all of the particles with that mass value would 

have that same velocity, the air-entry velocity would 

1 

be 12 km/sec, i.e., (6 2 + ll 2 ) 2 . But, as is likely, 
if V has a distribution with lower and higher 

values, the lower values are favored by tlie larger 
effective interaction cross section of the earth. 

The resulting air-entry velocity is then likely more 
nearly 11 than 12 km/sec. Then, when one takes 
11 km/sec for tlie average air-entry' velocity and 
combines it vectorially with that of a low-flying 
satellite (of about 7. 5 km/sec) , the average , 

velocity of impact is 13. 3 km/sec, i.e., (ll 2 + 7. 5 2 ) 2 . 
Then Dalton's [16] value for tlie meteoroid mass 
corresponding to tlie centroid of the Explorers XVI 
and XXIII puncture data should have been increased 
by a factor 10°' S7 ° instead of 10 0,283 . The resulting 
estimate of the mass of the just-puncturing meteoroids 
for tlie centroid of the Explorer data would be 
IQ-i.eo g, mos t probably a typical zodiacal light 
particle. 

It is difficult to judge for "B" whether or not 
tlie change in average air-entry velocity V from 

18.83 to 11 km/sec corresponds correctly with the 
correlation 0. 191 between mass m^ and velocity' 



V with terrestrial weighting, or whether or not 

CO 

the change in average geocentric velocity V from 

IS. 04 to G km/ sec corresponds with the correlation 
0. 188 between mass m^ and velocity with 

limar weighting, when mass is extrapolated down 
seven orders of magnitude. But the explanation is 
probably correct. This means that the Explorer 
point marked for sub-sample "B" in Figure 43 should 
be shifted 0. 283 units further to the right, bringing 
it into even closer agreement with the model extrapo- 
lated from the photographic meteor results. The 
Pegasus point would also be shifted somewhat to the 
right, but by a smaller increment due to the smaller 
difference in velocity from that of the photographic 
meteors. The point is sufficiently important that 
the satellite puncture data are compared in Figure 
44 with the model for the energy-cumulative influx 
extrapolated from the photographic meteor results 
in equation (52) . Again the agreement between the 
satellite puncture data and the model extrapolated 
from the meteor data is considerably closer with 
sub-sample "B" than with "A." Furthermore, for 
either "A" or "B, " the agreement with the satellite 
puncture results is considerably better when the 
extrapolation is with l'espect to energy in Figure 44 
than when it is with respect to mass in Figure 43. 

This tends to confirm the suspicion that the lack 
of closer agreement with the model extrapolated 
with respect to mass is due to the fact that the 
velocity distribution of tire meteor particles was 
erroneously presupposed for the satellite puncturing 
meteoroids in Figure 43. By shifting the satellite 
points in part (b) for "B" in Figure 43 to the right 
for the increment 0.2S3 in log m corresponding 
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to a change in mean impact velocity from 18. 83 to 
13.3 km/sec, the satellite points acquire the same 
configuration with respect to the mass model as they 
have in Figure 44 with respect to the energy model. 

Without the supporting evidence from the satellite 
puncture experiments, the extrapolation of the photo- 
graphic meteor models into tire mass and energy 
regions characteristic of the zodiacal light particles 
would have been quite eclectic. Nevertheless, one 
can consider the extrapolations well substantiated 
down to about 10 -6 g (which, according to Smith 
[3S] , is the mass of a typical individual grain in a 
dustball meteoroid) and a further extrapolation down 
to about 10 -7,6 g reasonably corporeal. Regardless 
of whether the large dustballs are derived from the 
smaller ones by accretion, or the small ones are 
derived from the larger ones by disruption, the 


mass-cumulative distribution should begin to "level 
off" toward the vicinity of the minuscule particles. 
The process of "leveling off" should be quite gradual 
due to the natural distribution of grain size. A 
suggestion of this effect is evident in Figure 44 in 
that the slope of the line joining the Pegasus and 
Explorer points is arithmetically less than that of 
the extrapolated model. The same trend was 
evident in the Explorer data which is represented in 
Figure 44 by its centroid. Further confirmation of 
this suspicion is given by the reported results of 
analysis of the zodiacal light and solar F-corona. 

In 1947 van de Hulst [98] found, by assuming that 
the log cumulative space density of zodiacal light 
particles is a linear function of log particle radius 
with slope - (p - 1) , that p has the value 2. 6. In 
1963 Ingham [72] reported the higher values 4 or 5 
for p derived by extrapolation from meteor data. 
Butin 1967 Blackwell, Dewhirst and Ingham [99] 
reported ”. . . it seems likely that the van de Hulst 
[98] density [concentration of particles] is the more 
representative of the optical determination of inter- 
planetary' space density. " Not enough is known 
about the relation between particle material density 
and mass over the range of zodiacal light particles 
to support an accurate formulation of the relation 
between the parameter p in the radius distribution 
and the parameter s in the mass distribution, 
equation (37) . But ignoring any dependence between 
radius and material density would give the value 3 
for (p - 1) / ( s - 1). Then, with the spatial weighting 
w in Figure 42, the values 2.42 and 2. 37 for the 

parameter s from "A" and "B, " respectively, 
would give the values 5.26 and 5. 11 for p from "A" 
and "B, " respectively'. The value for p from either 
"A" or "B" is somewhat high in comparison with the 
zodiacal light results, especially van de Hulst 1 s 
[98] ; but there is somewhat closer agreement from 
"B" than from "A. " 

Figure 43 illustrates models for the mass- 
cumulative terrestrial influx extrapolated by equation 
(50) to mass values below the range of meteor data. 

It is also of interest to assess the validity of a 
similar extrapolation of the results in equation (50) 
with lunar weighting Wj. By equation (50), at 10 _7 ’ c g 

in the models from sub-sample "A," the terrestrial 
influx exceeds the lunar influx by the factor 10~°' 7 ° 
or 5; and the corresponding factor for "B" is 
iO -0 ' 92 or S. As already noted, Opik's [97] value 
for the mean geocentric velocity' V of these small 

particles is 6 km/sec. If all of the particles had 
the same velocity', then the factor by which the 



terrestrial influx is increased by gravitational 

attraction would be 1 + (11/6) 2 or 4. 36. But one 

must expect that the particle velocities are distributed 

about their mean, with some values sufficiently 

small to increase the gravitational enhancement factor 

considerably. For example, any particles with a 

value 3 km/ sec for V would have a gravitational 
G 

enhancement factor of 14.44 but a particle with 
velocity higher than the mean by the same amount 
would have a gravitational enhancement factor of 
only 2.49. The average of the enhancement factors 
for 3 and 9 km/sec is 8.47, which is nearly twice 
as large as that for the G km/ sec particles. Thus, 
both values for the gravitational enhancement 
factor, 5 for the extrapolation from sub-sample 
"A" to 10 -7 ' 6 g and 8 for the extrapolation from "B, " 
are plausible; but the value 8 from "B" is more 
probable than the value 5 from "A." Anyway, 
nearly the same conclusions as for the mass- 
extrapolation of the terrestrial influx also apply in the 
extrapolation of equation (50) for the lunar influx; 
i.e. , one can consider the extrapolation fairly well 
substantiated down to about 10~ G g, with the further 
extrapolation down to about 10~ 7-6 g somewhat more 
diaphanous. 

The extrapolation of equation (52) for the 
energy-cumulative terrestrial influx is supported in 
Figure 44 by the satellite puncture results. The 
corresponding extrapolation of the energy-cumulative 
lunar influx in equation (52) is similarly supportable, 
but needs some explanation. By extrapolation of 
equation (52) , the lunar and terrestrial energ y- 
cumulative flux values are in the fixed ratios 10° -48 
or 3 and 10°' 32 or 2 from sub-samples "A" and "B, " 
respectively, independently of the extrapolation. It 
was already mentioned that meteoroids of 10 -7 ’ 6 g 
and average low-earth-orbit-satellite-impact 
velocity 13. 3 km/sec should be 5 to 8 times less 
abundant in the lunar influx and with a mean 
geocentric or lunar-orbiter-vehicle-impact velocity 
V about 6 km/sec. Because the particle energy' 

Li 

is reduced by a factor 5, and the influx is reduced 
by a factor 5 to 8, one might suppose the energy- 
cumulative influx to be reduced by a factor 25 to 40. 
But the low velocity particles, which have the high 
gravitational flux enhancement factors, have little 
effect on the cumulative distribution of energy. The 
punctures are made not so much by particles with 
average velocity and just-sufficient mass, but more 
by smaller particles with above-average velocity. 

What is significant is not the expected value of V , 

vjr 

but the expected value of V 2 ; 42 as in equations (44) 
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and (45) . This explanation and the extrapolation of 
the results from sub-sample "B" by equation (52) 
agree exactly with Gurtler and Grew's [100] results 
for the Lunar Orbiter puncture experiments. They 
(Gurtler and Grew [100]) reported that data collected 
between lunar altitudes 30 and 6200 km . . con- 
tinuously' for 17 months indicate (by 22 punctures) 
that the rate of penetration in the lunar environment 
is approximately half the rate in the near-earth 
environment as measured by detectors of the same 
type aboard Explorers XVI and XXIII. " The agree- 
ment with the corresponding extrapolated results 
from sub-sample "A" is somewhat poorer in that 
the "A" results indicate that the puncture flux 
reduction should have been by a factor 1/3 instead 
of 1/2. 

Comparisons to Other Reported Results 

As noted in 1966 by Nilsson [101] , for 
several years there has been a considerable volume 
of literature concerning the alleged existence of a 
concentration of dust in the vicinity of the earth, as 
inferred primarily from the satellite piezoelectric 
microphone measurement reported by Dubin and 
McCracken [102] and by Alexander, McCracken, 
Secretan, and Berg [103]. They had reported that 
the mass-cumulative impact rate, in particles per 
square meter per second for mass m in grams in 
the interval 10~ l ° ms 10 _G g, were 10 -17 ‘° m _1 ' 7 °. 
Nilsson's [101] paper in 1966 cast considerable 
doubt on the validity of those measurements by 
reporting satellite and laboratory evidence that the 
microphone crystals emit noise when subject to 
slowly varying temperatures. In 1968 Konstantinov, 
Bredov, and Mazets [104] reported experiments 
from the Soviet artificial earth satellite Kosmos-135 
which verified Nilsson's [101] contention. They 
reported using piezoelectric sensors, analogous to 
the others which had given the anomalously high 
readings, but specially developed and applied with 
careful precautions so as to protect the system from 
interferences of acoustic, thermal, or electrical 
origin. They reported good agreement with the 
puncture data from Explorers XVI and XXIII and 
Pegasus 1 and 2. Somewhat later, Koptev [105] 
wrote further about those experiments, which he 
said, although producing readings several thousand 
times smaller than anticipated, were in conformity 
with the American satellite puncture experiments 
with gas-filled pressure cells on Explorers XVI and 
XX11I and with electrical condensers on Pegasus. 

It has already been indicated that those puncture 



results agree well with the energy-cumulative 
distribution extrapolated by equation (52) from the 
photographic meteor results, the agreement being 
somewhat closer with sub-sample "B" than with 
"A. " 

In the present analysis die total flux (cumulative 
with respect to mass, momentum, or energy) was 
found by multiplying the sporadic flux by the factor 
2529/20G6 (or 1.22) by which the total exceeds the 
sporadics in McCrosky and Posen's [18] sample 
of photographic meteors. This may have constituted 
a slight over-correction, because the proportion of 
stream meteors which could have qualitied by the 
selection criteria for the sub-samples "A" and "B" 
might have been lower than the proportion of sporadic 
meteors which qualified. This aspect of the relative 
frequency of sporadic versus stream meteors was 
treated recently in an interesting analysis by Erickson 
[lOG] . He found that the total luminous flux of known 
streams is about 30 percent of the sporadic flux, 
that the total mass influx of known streams is of the 
order of one percent of that from sporadic meteors, 
and that the target puncture peak rates of the most 
active streams are no larger than the average 
puncture rate from sporadic meteoroids. 

The distribution of meteor-forming, vehicle 
puncturing, and zodiacal-iight-distributmg particles 
in interplanetary space may be quite different from 
the distribution of the bodies responsible for the 
craters apparent in the telescopic or Lunar Orbiter 
pictures of the lunar surlaoe and the Mariner IV 
pictures of Mars. According to an analysis by 
Anders and Arnold [ 107] , at least 75 percent of 
lunar craters down to a diameter of 20 km are 
caused by impact of bodies of asteroidal rather than 
of cometary origin, and the rate of crater formation 
on Mars by these asteroidal bodies is about 25 times 
higher than that on the moon. 

With respect to heliocentric distance r, by 
equation (33) the concentration in interplanetary 
space of the meteoroids (spatially weighted) from 
sub-samples ''A' 1 and "B" varies with r~°' 9, and 
r -1 ‘ respectively. Briggs [74] computed the 
spatial variation of particle concentration which 
would support the measured values of the zodiacal 
light under sleadv-state conditions of Poynting- 
Robertsor. retardation of debris continuously ensuing 
by the disintegration of comets, but exemplified by 
the appropriately weighted distribution of photographic 
meteor particles. Briggs [7-1] found that concentra- 
tion should vary with a power of r which would change 


from -0.97 through -1.09 to -1.27 as r changes 
from 0. 5 through 1. 0 to 1. 5 AU. His result for 
die exponent agrees well with the value -0. 97 from 
"A" in equation (33) for 1 AU. This result gives 
some support to the methods of the present analysis; 
but Briggs [74] used Hawkins and Soutlnvorth's [7G] 
random sample of 359 photographic meteors as an 
input with the mass values which were computed 
from "Model-A" rather than from "Model-B. " In 
19G7 papers by Soulhworth [108] and Divari [109] 
analyzed the space concentration of radio meteor 
particles. Divari [109] used the bivariate distribu- 
tion of orbit inclination and semi-major axis from 
7511 radio meteors by Kashcheev, Lebedenets, and 
Lagutin (see Divari [109]). statistically weighted 
to correct for the probability of earth encounter. 

He found that, for heliocentric distance r in the 
interval 0.5 s r £ 1.5 AU space concentration in 
the ecliptic plane varies nearly as r _f . Southworth 
[ 1 08] used data for 13 672 meteors observed from 
1962 to 1965 by the six-station radar network in 
Havana, Illinois. The meteors with negative 

celestial latitude of radiant B were discarded 
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and the remaining ones were carefully weighted for 
probability of earth encounter and other physical 
selective effects. The mass values in Soutlnvorth's 
[10S] sample ranged from about 10 -5 to 10 -3 g with 
optimum representation from 10 -3 - 5 to 10" 4 g. He 
plotted the logarithm of the relative space density 
versus the logarithm of the heliocentric distance 
r in the vicinity of the ecliptic by a series of 
connected line segments, evidently representing 
the 20 discrete stratifications made in r. Between 
Mercury and Earth the plot can be closely repre- 
sented by a line with slope -0. 77, and between 
Earth and the middle of the asteroid belt at 2. 7 AU 
by a line with slope -1. 5G. In addition to this 
"dog-leg” at 1 AU in the plot, there are also "dog- 
legs" seeming to hinge at the points for Mercury 
and Jupiter. These three appreciable bends in the 
plot at Mercury, Earth, and Jupiter may be spurious. 
The value -1. 56 for the slope between Earth and 
Jupiter agrees well with the value -1. 52 for sub- 
sample "B" in equation (33) ; but the mean (-1. 17) 
of the two values -i. 56 and -0.77 at the Earth agrees 
somewhat better with the value -0. 97 for "A" from 
equation (33). 

In an editorial reporting on the mission of the 
Mariner XV voyage to Mars in the 10 September 1965 
issue of "Science," Abelson [110] noted that the 
micrometeorile study shows that interplanetary 
dust was more abundant in the vicinity of Mars than 
near Earth. In the same copy of "Science" with 





Abelson's [110] editorial were a report by Anderson 
[111], describing the Mariner IV spacecraft, and 
a report by Alexander, McCracken, and Bohn [112] 
on the dust experiment. Impact measurements of 
cosmic dust were made on the Mariner IV space- 
craft throughout the nearly eight months of flight 
from Earth in November 1964 to Mars in 1965. 
Alexander, McCracken, and Bohn [112] reported 
an analysis of the 215 impacts of particles for which 
momenta were larger than 1. 96 x 10 -3 dyne-sec. 

The detector was described as an acoustical trans- 
ducer bonded to an aluminum sensor plate with a 
viewing solid angle of about n steradians. They 
reported an enhancement of the flux by a factor 5 
as the heliocentric distance r from the sun 
increased, and a change in the mass-cumulative 
impact flux from ~ m -0 * 55 near both planets to 
~ m -0-9 over the interval 1.25:= rs 1.36AU. 

But the reports did not mention any consideration 
of possibly changing orientation of the impact 
detector during the mission. Anderson [111] noted 
that the impact detector was on the body of the space- 
craft, with the microphone plate approximately 
perpendicular to the plane of the orbit when the 
axis of the guide-star Canopus was approximately 
perpendicular to the roll axis aligned with the 
heliocentric radial vector. But the reports did not 
mention the orientation of the detector with respect 
to tire Sun, nor whether or not the same guide-star 
was used throughout the flight. An analysis with 
respect to the detector plate orientation, similar 
to Dalton's [113] method for the interpretation of 
the orientation information in the puncture data 
from the Pegasus spacecraft, might account for 
the impact flux variations during the flight. More 
likely, however, the results are not valid if the 
sensor was of the type which Nilsson [101], 
Konstantinov, Bredov, and Mazets [104], and 
Koptev [105] have shown to give spurious signals. 

Notwithstanding the difference between "Model-A" 
and "Model-B" for photographic meteors, the bias 
and factors of physical selectivity in radio meteor 
data are less accurately known. The results are 
often quite disparate. For example, Figure 10 
illustrates that half of the meteoroids have radiants 
within • 26 degrees or ± 23 degrees of the ecliptic 
by sub-samples "A" and "B, " respectively; but 
Davies and Gill [60] noted in 1960 that half of the 
statistically weighted radio meteors have radiants 
within ± 15 degrees of the ecliptic. In 1963 Greenhow 
[ 114] discussed the peculiar difficulties of interpre- 
tation of radar meteor data in some detail and 
stressed the general lack of substance in prior 


results and the need for caution in future experiments 
and analysis. The most obvious difference between 
sporadic radio and photographic meteor data is in 
the velocity distributions. By the differences 
between the velocity distributions of sporadic 
meteors from photographic versus radio measure- 
ments, one does not imply that the velocity of an 
individual meteor would be assessed with such 
different results by the two methods, but that the 
samples are differently constituted by physical 
selectivity. This was illustrated in 1955 by 
Whipple [115] by tabulating radio and photographic 
measurements of the velocity of meteors in several 
distinct streams. A thoughtful paper by Furman 
[116], translated from the 1966 Russian publication, 
facilitates a realistic perspective for the difficulties 
remaining in the interpretation of radio meteor 
phenomena and for the uncertainty extant after 
careful analysis. Even more recently (in 1967) 

Fialko [96] has noted that so far an accurate 
dependence of meteor ionization probability p^ on 

velocity V has not been obtained (where ionization 

probability p^ means the probability that a single 

evaporated meteor atom will produce a free electron) . 
The hope of getting a more accurate determination 
of the velocity dependence of the luminous efficiency 
/ 3 through a determination of the velocity dependence 
of p and of the ratio fi/p does not seem to be very 

q q 

helpful so far. For example, in 1956 Hawkins [117] 
reported that P/p^ were ~ V~ 4-56 , and in 1964 

Hawkins [6] noted that a detailed calculation gives 

approximately p ~ V 3-4 . Solving for p between 

Q oo 

Hawkins [117, 6] 1956 and 1964 results would give 
approximately p ~ V -1-16 , in better agreement with 

"Model-B" than with "Model-A." But in both papers 
Hawkins [117, 6] gave the relation p ~ V as in 

"Model-A," and likely intended it to be understood 
as a more substantial result. On analysis of the 1956 
data on meteor echoes and magnitudes by Millmann 
and McKinley [118], Delcourt [119] deduced in 1963 
that the ionization coefficient jS (or ionization 

efficiency in the conversion of meteoroid kinetic 
energy into the freeing of electrons) is approximately 
related to p by P^/P ~ V^. But in 1968 Derbeneva 

[120] gave an analysis of experimental ionization 
and excitation cross sections for collisions of 
meteoroid atoms with nitrogen molecules and con- 
cluded that p /p is independent of V^. The onus. 
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evidently, is that Fialko's [96] notion in 1967. that 

so far an accurate dependence of p on V has not 

q «> 

been established, is correct. In 1967 Slattery and 

Friichlenichl (It’ll measured p for micron-size 

q 

iron particles injected at from 20 to 45 km/ sec 
into air targets, getting p ~ V 3 ' 12 , Their result 

is even more difficult to reconcile with ''Model-A" 
than was Hawkins' [6] ,1 ~ .[n comparing 

results from photographic and radio meteors, it 
should he noted that the photographic meteors span 
a transition region for radio meteors. The median 
mass for both sub-sampies " A ' ' and "B" is near 
the value 0. 1 g, which McIntosh [122j shows is the 
value below which, for a single non-fragmenting 
meteoroid, it is expected that log echo duration of 
radio meteors is a linear function of log' mass with 
slope (but not intercept) independent of velocity. 

The papers concerning the statistical dependence 
of meteoroid mass and velocity offer a clearer 
choice between "Model-A." with mass m. in sub- 

A 

sample "A" correlated only 0.020 with velocity V , 

and "Model-3, " with mass in sub-sample "B" 

correlated 0. 191 with V . Some of these papers, 

mentioned below, give a basis for a choice between 
"A. " with a low velocity distribution averaging 
16. 17 i 0. 17 km/ sec, and "B, " with a higher 
distribution averaging 18.83 ± 0.32 km/ sec. 

In 1955 Whipple and Hughes ! 123] noted that the 
average velocity of photographic meteors must be 
higher than the value 16. 5 km/ sec which they 
adopted for the air-entry velocity of meteorites. 
This condition is satisfied by ”3, " with average 
velocity iS. 83 km/sec, but not by "A" with average 
velocity 16.17 km/sec. 

In two papers presented in 1957, Whipple [124, 
125] noted that "... a velocity of 28 km/sec is 
average for photographic meteors. Undoubtedly 
the velocity falls off for smaller meteoroids as we 
deal more and more with particles whose orbital 
eccentricities and dimensions have been reduced 
by physical effects. . . A mean value of 15 km/sec 
has been arbitrarily chosen for the smallest 
meteoroids and an arbitrary gradation of velocity 
with magnitude adopted. " The value 28 km/sec 
is near the geometric mean velocity for a random 
sample of photographic meteors, whereas the mean 
values 16. 17 and 18.83 km/sec lor "A" and "B. " 


respectively, are for random samples of meteoroids 
in the terrestrial influx. Those results favor "B" 
over "A" both in the higher mean value (18. 83 for 
”B" versus 16. 17 km/sec for "A") and in the higher 
correlation (0. 191 for "B" versus 0.020 for "A") 
necessary to account for the change in mean velocity 
from 2S to 15 km/sec as mass decreases. 

In 1964 Hawkins, Lindbiad, and Southworth [126] 
described a controlled experiment in which the veloc- 
ity distribution of radio meteors was measured at 
several different limiting magnitudes, with minimi- 
zation of various selection effects. Mean aii'-entry 
velocity of all meteors brighter than visual magnitude 
+ 9 was found to be only 87 percent of that for all 
meteors brighter than visual magnitude -5-6; i.e. , 

33.3 versus 3S. 3 km/sec. The velocity measure- 
ments were considered preliminary. They are 
subject to the continuing disparity between photo- 
graphic and radio assessments of the statistical 
distribution of meteor velocity. Balter [127] had 
considered the average velocity decrement in that 
experiment to be only 3. 9 instead of 5. 0 km/sec, 
but concluded that the data still indicate a decrease 
in average meteor velocity with fainter meteors. 

The effect was attributed to the difference in orbits 
within the various meteor populations, in agreement 
with results attributed by McKinley [5] to Opik. 

Also Hawkins [128] suggested that Eshleman and 
Gallagher [ 129], who found in 1962 that the average 
velocity of sporadic meteors decreased by 2 km/sec 
between magnitudes +7 and +12, would probably have 
noted a larger decrement in mean velocity’ if they 
could have included values below 35 lon/sec. 

Although the disparity in the velocity distribution 
of radio versus photographic meteors makes one 
apprehensive about the integrity of these results to 
the extent that they might be overlooked except that 
they do agree with other results, they support 
choosing "B" for its significant correlation 0. 191 
between mass and velocity V^. 

In 1964 Opik [97] noted that the average geo- 
centric velocity V for dustballs and zodiacal 

light particles can be assumed equal to 18 and 6 
km/ sec. From Figures 15 and 16 with lunar weight- 
ing, the values for the mean velocity' and correlation 
between mass and velocity in the lunar influx are 
13.66 km/sec and 0.024, respectively', for "A" and 
IS. 04 km/sec and 0. 18S, respectively, for "B. " 
Opik’s [97] values give a clear support for choosing 
the "B" results. 
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In summary of his extensive literature survey 
and study of meteors in 1966, Vedder [ 8 ] noted that 
. . the fainter meteors seem to have lower 
velocities, possibly a results of the Poynting- 
Robertson effect. . . " and that the calculation of the 
masses and flux is complicated by the effects of 
fragmentation and radiative cooling. This supports 
choosing "B" for its significant correlation between 

mass in and velocity V . 

B 

In 1967 Lebedinets and Kashcheev [130] reported 
statistical results for radar measurements of 
velocities and radiants and corresponding orbit 
computations for 12 500 meteors brighter than 
magnitude +7. The results, which they said were 
corrected for physical selectivity’ with respect to 
velocity and radiant position, were reported in com- 
parison with corresponding results by McCrosky 
and Posen [18] for Super-Schmidt photographic 
meteors brighter than magnitude +3. They 
(Lebedinets and Kashcheev [130]) found that the 
radio meteor sample, with a threshold four mag- 
nitudes fainter (or about 1/60 lower mass threshold) 
than the photographic meteor sample, gave sub- 
stantially different distributions of orbital elements. 

It was found that by considering meteoroids of 
smaller and smaller mass there is a systematic 
decrease in the mean dimensions of the orbits. 

The mean perihelion distance q for the radio meteor 
particles was found to be approximately half as 
great as for the photographic meteors. No comment 
was offered that these results imply a lower mean 
air-entry velocity for the smaller particles. 

Velocity (e. g. , as is shown by Opik [70]) 

depends also on eccentricity e and the inclination 
i of the orbit to the ecliptic, and the authors showed 
that the distribution of i, too, changes with particle 
mass. But the results lead one to suspect that the 
smaller particles may have a lower mean air-entry 
velocity V . This interpretation is substantiated 

by Kresak's [20] carefully weighted analysis also 
of the data of McCrosky and Posen [18] . By these 
results, by Lebedinets and Kascheev [130] and 
Kresak [20], one expects some positive correlation 
between mass and velocity, as in "B, " rattier than 
"A. " Further basis for that expectation is given 
by Ceplecha's [45] 1968 analysis of McCrosky and 
Posen's [18] photographic meteor sample. He 
found that the proportions of meteors in four 
distinct groups vary with brightness and velocity, 
and that the characteristics of the orbits vary from 
group to group. McCrosky and Posen [18] noted 


also that the orbital characteristics of photographic 
meteor particles vary considerably with meteoroid 
size. 

A number of published works treating the 
distribution of meteor particle velocity and its 
dependence on mass have been mentioned. Their 
interest with respect to a choice between "Model-A" 
with sub-sample "A" and "Model-B" with "B" is 
quite indirect and not fully satisfying. Beginning in 
1956 with Levin's [3] comprehensive treatise, a 
number of works have appeared which offer studied 
opinion, creative insight, unusual perspective, or 
accidental involvement in the ambit of meteor 
luminous efficiency. Some of these works will be 
mentioned further. The "Model-A" tends to be 
regarded as a truism after three decades from its 
first use by Whipple [25] in the determinations of 
the density of the upper atmosphere with an accuracy 
which has been supported subsequently fairly well 
from the behavior of orbiting satellites and re-entry 
bodies. But, as Levin [3] noted specifically’ and 
as Lovell [84], Babadzhanov [i31], and Katasev 

[132] formulated also in mathematical detail, the 
calculated atmospheric density p is inversely 

proportional to the (1/3) -power of the meteor 
luminous efficiency /3 (and a function of other 
parameters) . Therefore, the low accuracy of /l 
is only slightly reflected in p - but the meteoroid 

mass m ~ has no gain in accuracy. Both 

Delcourt [119] in .1963 and especially Kovshun [133] 
in 1966 considered "Model-A” to be only an approxi- 
mation and expedient to the detailed and thorough 
theory of meteor luminosity which in die present 
analysis is approximated by "Model-B. " Kovshun 

[133] noted that for direct determination of d in each 
specific case the formulation approximated by "Model- 
B" ". . . is extremely inconvenient . . . For this reason, 
meteor mass is computed to this day from the 
simplified theory of meteor luminosity [i.e. , "Model- 
A 1 '] which does not allow for the dependence of on 
the mass, structure, and composition of meteoroids. " 

In 1967 Jacchia, Verniani, and Briggs [56] 
reported a statistical analysis of the in-atmosphere 
characteristics of 413 precisely reduced photographic 
meteors. They compared theoretical and empirical 
relations presupposing the "Model-A" luminous 
efficiency ft = TnV and corresponding mass value 

m^. It is of interest to note how some of the 

empirical results would seem to be brought into 
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closer agreement with the theory if the "Model-B" 
luminous efficiency P and corresponding mass 

-D 

nig were replacing P ^ and , respectively. 

Theoretically the meteor beginning height for 
meteoroids of a given mass should depend on the 
constant value of (?pV 5 , but the empirical 

CO 

criterion is pV 3-5 , where p is the atmospheric 

CO 

density. The velocity exponent is unduly reduced 
by using p^ instead of p for p. The maximum 

photographic luminous intensity I ^ was found to 

be ~ V 3 ’ 5 (i. e. , 8. 75/2. 5) instead of the theoretical 

°o 

value ~ /3V^. Thus, the velocity dependence of p 

would be more nearly ~ V 0 - 5 than in "Model-A" 
witli ~ V . On tlie average this result would agree 

CO 

somewhat more closely with "Model-B. " The 
empirical beginning height criterion pV^ s found 

by Jacchia, Verniani, and Briggs [56] agrees quite 
closely with Levin's [134] velocity dependence of the 
ratio of apparent radiant density of visual meteors 
and true radiant density, ~ V^ 6 . 

In one of the Smithsonian contributions to 
astrophysics in 1965, Verniani [26] equated the 
dynamic and photometric mass values in the data 
subsequently reported as "An Analysis of the 
Atmospheric Trajectories of 413 Precisely Reduced 
Photographic Meteors," by Jacchia, Verniani, and 
Briggs [79], In his analysis presupposing an 
exponential velocity dependence of luminous 
efficiency p, Verniani [26] concluded that he had 
demonstrated the verity of p ~ V per se in 

"Model-A, " independently of its original basis and 
limitations (bright meteors only) in accordance with 
the model adapted by Whipple [25] in 1938 from 
Opik's [19] 1933 "Atomic Collisions and Radiation 
of Meteors. " Verniani [26] implies that Jacchia's 
[90] finding (confirmed by Jacchia, Verniani, and 
Briggs [56]) , that meteor color index at constant 
magnitude is independent of velocity, should support 
the assumption that luminous efficiency is propor- 
tional to an exponential function of velocity. He 
showed that the independence of velocity and color 
index, which is proportional to the logarithm of the 
ratio of the photographic and visual luminous 
efficiencies, establishes that the photographic and 
visual luminous efficiencies can differ only in their 
coefficients (the luminous efficiency factors) and 
not in their velocity exponents. Actually Jacchia's 
[90] discovery has a broader implication, i.e. , 


that the photographic and visual luminous efficiencies 
(for meteors of nearly the same magnitude) are 
proportional to the same function of velocity, 
whether exponential or not. It does not appear that 
Verniani's [26] discussion on this point has any 
bearing on the question of whether or not luminous 
efficiency might be adequately approximated as an 
exponential function of velocity. By equating the 
dynamic and photometric mass formulations, 

Verniani [26] got an explicit expression for P/ p 2 ^, 

the ratio of the luminous efficiency and the square 
of the meteoroid material density. He attributed 
all of the variation in j S/p 3 to the variation of p 

with velocity V and discarded the 10 percent of 

the sample with the largest deviations of P/p 2 m 

from its mean. In 1967 Verniani [135] repeated 
his analysis of /3/p 2 ^ with the same sample he had 

used in 1965 [26]; but this time he discarded 22 

percent of the 247 sporadic meteors which he had 

not already discarded before (leaving only 189 for 

his new analysis) . This was done by presupposing 

p ~ V and discarding the meteors with p/p 2 
» m 

differing most from its consequent mean. Thus, 
the result which Verniani got, P ~ V*‘ ^ ^ in 

OO 

1965 [26] and P V^ - ° ~ °' 1 in 1967 [135] , is 
probably mostly a reflection of his two presupposi- 
tions that P ~ V° and that n is unity. 

Verniani's works [26, 135] have been widely 
referenced as evidence of the essential verity of 
"Model-A" with respect to the velocity variation of 
meteor luminous efficiency p = t 0 V^ and as an 

authoritative refinement of the value of the constant 
luminous efficiency factor t 0 ; i. e. , 10 _0 ’ os or 81 
percent of the value of r 0 in "Model-A. " For 
example, Verniani [26] stated that in 1963 Whipple 
[136] used Verniani's value -18. 91 for log r 0 p 

(subtract 0. 72 to convert to visual units, add 9. 72 
to convert from Jacchia's magnitude units to CGS 
units, and add 5. 00 to convert from cm/sec to 
km/sec velocity units) in the construction of the 
mass-cumulative influx model for meteors, based 
on Hawkins and Upton's [81] magnitude cumulative 
influx model. The value -18. 91 was based on an 
interpretation of the Trailblazer I artificial meteor 
rocket experiment by McCrosky and Soberman [137], 
Verniani [26] stated that the reliability of that value 
based on Trailblazer I has to be .considered higher 



than that obtained with three "asteroidal" meteors 
by Cook, Jacchia, and McCrosky [138] in 1963. The 
value from Trailblazer I was 1(T 0 - 42 , or 38 percent 
of the geometric mean of those by the three natural 
meteors studied by Cook, Jacchia and McCrosky 
[138], Whipple adopted the geometric mean of 
those two separately reported sets of values, 
getting a visual luminous efficiency factor 31 
percent of the value in "Model-A. " The value 
advocated by Verniani [135] in 1967 for tq is 81 
percent of that in "Model-A. " But in three recent 
papers Ceplecha [44, 43, 45] has spoken out in- 
creasingly more critically of Verniani's [26, 135] 
presupposition that a formula for /Vp 2 ^ could be 

derived by equating the expressions for dynamic and 
photometric mass values. In 1956 Levin [3] had 
criticized Whipple's earlier practice of equating 
the dynamic and photometric mass values to get 
inordinately low values for the meteoroid material 
density p . Levin [3] offered the explanation lhat 

". . . the minuteness of the dynamic masses of 
meteoroids is a result of their partial fragmentation, 
not of the minuteness of their densities. In the 
case of fragmentation, the mass determined by 
deceleration is the mean mass of the individual 
fragments, while the photometric data give their 
lotal mass. " 

In 1967 Allen and Baldwin [46] reported obser- 
vations, in ground laboratory facilities, showing 
that rock will froth during ablation, and that the 
froth may be sloughed as a result of aerodynamic 
pressure forces. Their results offer an alternative 
to Lhe very low material densities usually reported 
for meteors; but, of more pertinent relation to the 
present analysis, their experiments also checked 
the velocity dependence of luminous efficiency. They 
found, when luminous efficiency ft is considered tc 
be the product of a luminous efficiency factor tq 
and velocity V , that their experimental results 

showed that r 0 must decrease with increasing V . 

This result is contrary to "Model-A," which has r 0 
constant, but die variation of r 0 with V is in the 

sense that would tend toward agreement with "Model- 
B, " which on the average has a weaker velocity 
dependence than "Model-A," i. e. , the results by 
Allen and Baldwin [46] indicate that the velocity 
dependence of luminous efficiency is not as strong 
as in "Model-A. " A later study by Baldwin and 
Allen 1 139] confirmed Uiat photographic meteor data 
are consistent with their compact stone particle 
hypothesis, and that luminous efficiency ft depends 


less strongly on velocity V than that indicated for 

CO 

"Model-A." Two other NASA -sponsored laboratory 
measurements of the luminous efficiency of hyper- 
veiocity particles are pertinent to a choice between 
"Model-A" and "Model-B. " hi 1965 the Utah 
Research and Development Company [140] reported 
on their contracted research with sub-micron hyper- 
velocity spray particles generated by the hypervelocity 
impact of aluminum and steel balls into similar 
targets. With successful measurements under con 
trolled conditions, the spray particles with radii 
between 0. 005 and 0. l micron ar.o speeds between 
3 and 18.4 km/see, but with most of the data between 
6 and 16 km/ see, the luminous efficiency proved to 
be independent of speed. This result is contrary 
to "Model-A" but would agree considerably better 
with "Model-B" on the average. In 1968 
Friichtenicht, Slattery, and Tagliaferri [141] 
reported luminous efficiency measurements in the 

spectral band 3400 to 6300 A for sub-micron diameter 
iron particles injected from an electrostatic particle 
accelerator into air targets. The luminous efficiency 
ft was reported to be nearly constant over the velocity 
interval from 20 to 40 km/see, They stated further 
that ". . . within the limits of the experimental errors, 
ft actually appears to decrease slightly with increas- 
ing velocity. " The results of Friichtenich, Slattery, 
and Tagliaferri [141] also tend to be compatible with 
"Model-B," but not with "Model-A." 

In 1962 Levin 1142] noted thaf the presumption 
of meteor luminous efficiency ft proportional to 
velocity V independently of meteoroid mass m 

would have special implications for fragmenting 
meteors. For a meteoroid with given mass and 
velocity, the integrated luminous intensity would be 
independent of the meteor path-length L. The 
increase of the mean light intensity would have to be 
directly proportional to the decrease of the path- 
length. Detailed implications of that result, in 
effecting the parameter s in the mass distribution 
function, equation (37). were given in 1965 by 
Kruchinenko [143j. But in the same issue of the 
journal with Kruchinenko's ] 143] paper, Kramer 
[144] gave data showing that the theoretical 
predictions presupposing meteor luminous efficiency 
independent or mass, as in "Model-A, " do not con- 
form to reality in that the integrated radiation from 
a fragmenting metov -‘old is always less than that 
which would be inferred by Levin's [142] rule. The 
"Model-B" agrees with Kramer's [144] results 
better than "Model-A” in that luminous efficiency 
fty^ is a monotcnically increasing function of particle 

mass m over a -.vide range of values of m_. 

.D Id 



By the results of the present analysis, compari- 
son of consequences, and sampling of seemingly 
relevant material coming to the author's attention, 
the undeniably oversimplified "Model-A" must be 
rejected in preference for "Model-B" on grounds of 
relative plausibility of consequent implications. The 
undobtedly extreme complexity of the meteor 
luminous efficiency is substantiated by recent papers 
on meteor spectra, e. g. , by Hoffman and Longmire 
[145] and by Smirnov [146]. In context with such 
investigations, it becomes clear that one must not 
interpret the present analysis as any proof of the 
validity of Opik's [2] physical theory of meteors 
(approximated in the present analysis by "Model-B") , 
but only recognize that the inherent consequences 
led to no discernible contradiction or disagreement 
with other established fact so far as the analysis was 
pursued. On the other hand, the author feels that 
the case against any formulation of luminous 
efficiency proportional to velocity independently of 
mass, such as "Model-A," has been presented so 
sufficiently that it could be mollified only by 
legerdemain. 


m A for photometric mass. Therefore, the same 

value -2. 25 for the coefficient k 2 would have been 

found if mass m in equation (36) had been replaced 
A 

by /3 m the model-invariant product of luminous 
B B 

efficiency and photometric mass. This follows 

because m„ is the same product as fl. m. , and 
B B A A 

|8 is independent of m^. But p ^ is not independent 

of m . The definition of magnitude M, as introduced 
B 

in 1856 by Pogson (see Jones [147]) for stars, is 

sufficient for the proof 3 that mass m in the right 

B 

side of equation (37) can be replaced by /? m ■ i. e. , 

if B 


N 

>m 


(m) = N (1) 
>m 


,<l + kj) ( 


k 2 log x) 


(53) 


where k 2 is a constant in the exponential approxima- 
tion of p as a function of m , 

B B 


Now that "Model-A" has been eliminated in favor 
of "Model-B, " it is of interest to return to equations 
(36) and (37) and reconsider the role of the param- 
eter k 2 with the value -2.25 which Jacehia, Verniani, 
and Briggs [56J found by least-squares statistical 
analysis of photographic meteor data. In equation 
(36) the values for mass m used by Jacchia, 
Verniani, and Briggs [56] were the "Model-A" values 


k 2 

/S B ~m B 2 . (54) 

A check on the relation in equation (54) was made 
with a constant value for V of 19 km/sec, near 

OO 

the mean value 18. 83 ± 0. 32 km/sec for the terres- 
trial influx (Fig. 15) . A plot of log (j3 m ) versus 

B B 


3. 2.512 


-M 

" W 0 * 


••• M = (- log 2. 512) [log I M (M) - log 1 M (0) ] 

= kg + kj log + k 2 log (fi m) + k 3 log cos 
••• I ftI (M) ~ (p m) (M) 2 ' @ const. V^, etc. 


M 0m (P m ) ~ k 2 lo & 


(p m) M (M) 


F M [ M /3m (/i,U) ] = V° } X 


V m) 


W CO 

/ f m (m)dm= / f M [Mtfm)] dM tfm) ~ [(Pm) M (M)] 


k 2 log x 


/3 m 

( 1 + k 2 ) ( k 2 log x) 


where /3 ~ m 


k 2 
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log m over the interval 10 i m £ 10 g showed 

no graphically discernible deviation from linearity. 
The value of the exponent k 2 in equation ( 54) was 
found to be 0. 069 when V is 19 km/sec. This 

oo 

means that the value of the mass exponent - (s - 1) 
deduced from Hawkins and Upton's [81] analysis of 
meteor data is changed from -1. 21 to -1. 29. This 
value -1. 29 compares very favorably with the value 
-1. 31 found directly from the cumulative distribution 
of mass values in Figure 40. Although the value of 
kj in equation (54) was evaluated only for 19 km/sec, 
it is probable that the resulting value 0. 069 for k 2 
is somewhat lower than the effective mean value 
averaged over the range of velocity V . This is 

OO 

evident from Figure 4, which shows that the mass 
dependence of the luminous efficiency y3 is stronger 

for the higher values of velocity V . It is not un- 
reasonable to expect that the effective value of lc 2 
might be as high as 0. 10. If so, then, by equation 
( 53) , the entire discrepancy between the theoretical 
value -5/2 for k 2 and the empirical value -2.25 found 
by Jacchia, Verniani, and Briggs [56] would be 
accounted for through the factor (1 +1c 2 ) due to the 
dependence of meteor luminous efficiency p on 

meteoroid mass m n from "Model-B. " But the 

author's opinion is that the correlation 0. 191 ± 0.037 
between mass and velocity may play some 

part in the role of k 2 . Therefore, the value -1.21 
for the exponent - (s - 1) in the energy- cumulative 

distribution of meteoroids in the terrestrial influx 
is retained as in equation (52) and Figure 44b. 


CONCLUSIONS 


The greater merit of Opik's [2] 1958 physical 
theory of meteors, approximated in the present 
analysis by "Model-B" [equation (13) and Fig, 4], 
relative to Opik's [19] 1933 provisional theory 
[approximated by "Model-A," equation (11)], or 
any theory for luminous efficiency p proportional 
to velocity V independently of photometric mass 

OO 

m, is demonstrated by the present analysis. 


The meteoroid natural environment for the 

earth is, as expressed by F , F ,, , and 

>m ’ >m_V 
B B » 

F , , r2 , the mean cumulative influx per square 

> m_ r 

- B «> 

meter per second with respect to mass, momentum, 
and kinetic energy, respectively. 


log F = -1. 31 log m B - 14.77 (Terrestrial) 
>mB (50) 

l0gF >m V = -1-28 log m B V n 
B °° 

-13.07 (Terrestrial) (51) 

logF >A m V 2 = - 1 -2i i°g (im B V^) 
i B 

-11.72 (Terrestrial) , (52) 

where mass m_ is in grams and velocity V is in 

, / B 00 

km/ sec. 

The corresponding meteoroid natural environment 
for the moon is: 


log F >m = -i.21 log m B - 14.93 (Lunar) (50) 
B 


log F = -1. 28 log m V - 13. 38 (Lunar) 

>111--. V B Cj 

B G (51) 

l0gF >^m V 2 =- 1 - 21l0g ^ m B V G> 

B Cj 


12.04 (Lunar) 


(52) 


The meteoroid natural environment for an 
interplanetary vehicle at heliocentric distance r in 
AU differs from that for the moon only in that the 
particle mass-cumulative impact flux varies with 
r -2 ’ 05 , by equation (32) , and the mean impact 
velocity varies with r -0-50 . The concentration of 
particles in space varies as r _1 , equation (33). 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812 
September 3, 1969 124-091400 
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LUMINOUS EFFICIENCY OF 



FIGURE 1. ORIGINAL BASIS FOR THE HARVARD LUMINOUS EFFICIENCY OF METEORS 
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FIGURE 2. RELATION BETWEEN THE OPIK |2, 28] (1958, 1963) LUMINOUS EFFICIENCIES 
OF COMPACT IRON AND STONY DUSTBALL METEOROIDS 
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FIGURE 3 . RELATION BETWEEN THE OPIK [2, 28] (J 
OF COMPACT STONY AND STONE DUS 
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FIGURE 5. "MODE L-C" APPROXIMATING OPIK'S [2, 2S] (195S. 1963; LUMINOUS EFFICIENCY 
@ c FOR COMPACT STONE METEOROIDS £ lo 2 g 
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FIGURE 6a. MASS VALUES BY "MODEL-A" VERSUS "MODEL-B," SUB-SAMPLE "A" 
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FIGURE 6b. MASS VALUES BY "MODEL-A" VERSUS "MODEL-B, " SUB-SAMPLE "B" 
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FIGURE 7a. MASS VALUES VERSUS AIR-ENTRY VELOCITY, SUB-SAMPLE "A" 
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FIGURE 7b. MASS VALUES VERSUS AIR-ENTRY VELOCITY, SUB-SAMPLE "B" 
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CELESTIAL LATITUDE OF RADIANT, £ 


FIGURE 8a. ZENITH VERSUS RADIANT DEVIATIONS FROM THE ECLIPTIC, SUB- SAMPLE "A" 
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CELESTIAL LATITUDE OF ZENITH, /j 



CELESTIAL LATITUDE OF RADIANT, 


FIGURE Sb. ZENITH VERSUS RADIANT DEVIATIONS FROM THE ECLIPTIC, SUB-SAMPLE "B" 
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CUMULATIVE DISTRIBUTION WITH WEIGHTING, 



FIGURE 10a. DISTRIBUTION OF NEGATIVE VERSUS POSITIVE CELESTIAL LATITUDE 
WEIGHTED FOR SURVEILLANCE AREA AND APPARENT FRACTION OF CIRCLE 
OF CELESTIAL LATITUDE, SUB-SAMPLE "A" 
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CUMULATIVE DISTRIBUTION WITH WEIGHTING 



FIGURE 10b. DISTRIBUTION OF NEGATIVE VERSUS POSITIVE CELESTIAL LATITUDE 
WEIGHTED FOR SURVEILLANCE AREA AND APPARENT FRACTION OF CIRCLE 
OF CELESTIAL LATITUDE, SUB-SAMPLE "B" 
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CUMULATIVE DISTRIBUTION WITH WEIGHTING 









CUMULATIVE DISTRIBUTION WITH WEIGHTING 
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CELESTIAL LATITUDE WEIGHTED FOR 
SUB-SAMPLE "B" 








Weighting 

Uniform 

Terrestrial 

Lunar 

Spatial 



0. 267 

0. 298 


Mean Mass 


0. 207 

0. 225 

■EH 

Standard Deviation 

■■ 

0.564 

0. 583 

mam 

of Mass 


0.362 

0. 366 

ymMm 

Other Parameters 

Correlation Between Mass and Other Parameters, Times 10 3 

\r 

221 

191 

232 

167 

V 

CO 

212 

ISO 

227 

184 

u 

219 

188 

228 

169 

V G 

211 

177 

225 

1S6 


202 

186 

205 

222 

6 

210 

1S8 

222 

247 



,U A 

S05 

856 

776 

S21 

833 

875 

807 

823 

h b 

277 

249 

283 

262 
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248 

291 

285 


NOTE: Upper and lower entries in each table-block involve "Model-B” and ''Model-C" mass values, 
respectively. 

FIGURE 14. STATISTICAL ROLE OF MATERIAL CONSISTENCY, PRESUPPOSING STONE DUSTBALL 
VERSUS COMPACT STONE PARTICLES, BY MODELS FOR OPIK'S [2, 28] (1958, 1963; 

RESULTS WITH SUB-SAMPLE "B" 
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FIGURE 16. CORRELATION ARRAY - BLOCKS WITH TERRESTRIAL WEIGHTING (ABOVE THE DIAGONAL) 

AND WITH UNIFORM WEIGHTING (BELOW) 
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FIGURE 17. CORRELATION ARRAY - BLOCKS WITH LUNAR WEIGHTING (ABOVE THE DIAGONAL) 

AND WITH SPATIAL WEIGHTING (BELOW) 
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FIGURE 18b. APEX-RADIANT LONGITUDE DISTRIBUTION, TERRESTRIAL WEIGHTING, 

SUB-SAMPLE "B" 
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FIGURE 27b. PERIHELION PROBABILITY DENSITY FUNCTION. SUB-SAMPLE "B" WITH LUNAR 
WEIGHTING w.. AND SPATIAL WEIGHTING w 
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FIGURE 35b. DISTRIBUTION OF ARITHMETIC CELESTIAL LATITUDE OF RADIANT, 
LUNAR WEIGHTING, SUB-SAMPLE "B" 
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FIGURE 36b. DISTRIBUTION OF SINE INCLINATION OF ORBIT TO ECLIPTIC, SPATIAL 

WEIGHTING, SUB-SAMPLE "B" 
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FIGURE 39b. DISTRIBUTION OF AIR-ENTRY VELOCITY FOR TWO MASS REGIMES, 
TERRESTRIAL WEIGHTING, SUB-SAMPLE "B" 
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FIGURE 44b. ENERGY-CUMULATIVE INFLUX OF SATELLITE-PUNCTURING METEOROIDS VERSUS 
EXTRAPOLATED MODEL FROM PHOTOGRAPHIC METEORS, SUB-SAMPLE "B" 
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